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I.  INTRODUCTION 


(c)  During  November  1975,  the  Moored  Surveillance  System  (MSS) 

Field  Validation  Test  (FVT)  was  conducted  under  the  sponsorship  of 
the  MSS  Project  Office  (PME  124-30)  of  Naval  Electronic  Systems  Command. 
Applied  Research  Laboratories,  The  University  of  Texas  at  Austin  (ARL:UT) , 
participated  in  the  processing  and  analysis  of  he  acquired  data.  This 
report  contains  the  results  of  the  work  performed  under  Contract 
N00039-77-C-0003.  The  purpose  of  this  report  is  to  provide  measurements 
and  analyses  of  candidate  sensor  performance  based  on  data  from  the 
MSS-FVT, 


(U)  This  report  is  partitioned  into  four  volumes.  Volume  I describes 
the  data  collection  system  and  the  measurement  system  used  to  obtain 
these  results.  Volume  II  contains  the  detailed  data  products  obtained 
with  standard  frequency  resolution  processing.  Volume  III  contains 
the  detailed  data  products  obtained  with  vernier  frequency  resolution 
processing.  Volume  IV  contains  the  background  information,  summary 
data  products,  and  analysis. 

(U)  Since  the  MSS-FVT  was  completed,  the  name  of  the  MSS  program  has 
been  changed  to  Rapidly  Deployable  Surveillance  System  (RDSS).  Not 
only  the  name  but  also  the  scope  and  requirements  of  the  program  have 
undergone  major  revisions.  To  avoid  ambiguity,  the  term  MSS  will 
be  used  throughout  this  report;  however,  the  issues  addressed  herein 
are  those  of  current  interest  to  RDSS,  as  specified  by  the  project 
office  at  the  time  that  this  report  was  written. 


(U)  Volume  I is  intended  to  provide  a description  of  the  complete  data 
collection  and  measurement  systee  used  to  obtain  the  results  contained 
in  Volumes  II  through  IV.  The  contents  of  this  volume  are  divided  into 
seven  sections.  Section  II  describes  the  sensors  whose  performance  is 
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(U)  evaluated  by  this  report.  Section  III  briefly  describes  the  data 

collection  and  preprocessing  systems  that  input  to  the  ARL:UT  measurement 
system  (Fig.  1-1).  The  fourth  through  sixth  sections  describe  the  three 
major  components  of  the  ARL:UT  system  (Fig.  1-2):  A/D  Conversion  and 

Data  Reduction,  Ambient  Sound  Field  Measurement  Techniques,  and  cv 
Measurement  Techniques.  The  final  section  contains  recommendations 
pertaining  to  future  performance  analyses. 

(C)  The  sensor  performance  measurements  and  analyses  contained  in  this 
report  do  not  address  the  entire  MSS  processing  schetie.  This  effort 
vas  designed  to  directly  address  the  following  issues. 

• DIFAR  and  differenced  DIFAR  beamformer  performance  improvements 
over  an  omnidirectional  sensor. 

• Signal,  noise,  and  array  gain  measurements  for  an  MSS  candidate 
deep  moored  sensor. 

• Bearing  accuracy/performance  achievable  with  the  ACODAC  system/MSS 
candidate  acoustic  array. 

• Clutter  related  (line  loading)  statistics  for  a deep  moored  acoustic 
array. 

Issues  involving  signature  formation,  classification,  interbuoy  target 
formation,  localization,  field  management,  or  field  performance  are 
not  directly  addressed. 

(C)  At  the  time  that  this  work  was  Initiated,  the  ARL:UT  measurement 
system  did  not  contain  many  of  the  required  features  (i.e.,  DIFAR 
beamforming  and  clutter  measurements).  The  sponsor  provided  APX:UT 
with  an  outline  of  the  desired  measurement  system  characteristics  and 
a set  of  suggested  algorithms  (Ref.  1).  Much  of  the  time  and  resources 
expended  on  this  effort  were  used  to  implement  these  features.  It  was 
then  found  that  the  specified  line  formation  and  line  tracking  algorithms 
did  not  perform  well  on  these  data.  These  algorithms  were  subsequently 
revised  by  ARL:UT  because  a reasonable  performance  level  was  required 
for  the  execution  of  this  work. 
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FIGURE  l-l 

PREPROCESSING  DATA  PATH 
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FIGURE  1-2 

ARL.-UT  MEASUREMENT  SYSTEM  OVERVIEW 
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There  are  currently  In  existence  several  different  systems  which 
have  been  used  to  measure  sensor  performance.  Each  of  these  systems 
was  designed  with  different  constraints  and  goals,  and  their  designers 
frequv  tly  choi.e  to  use  significantly  different  algoiithms.  This 
report  is  intended  to  furnish  sensor  performance  measurements  which 
can  be  compared  with  those  from  other  measurement  systems.  To  help 
facilitate  such  comparisons,  the  ARL:UT  algorithms  end  their  associated 
assumptions  and  models  are  herein  described.  Since  the  ARL:UT  system 
functions  in  a laboratory  environment  without  the  usual  realtime  or 
tactical  constraints,  most  of  the  timing,  hardware,  and  formatting 
details  have  been  deleted. 
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II.  SENSOR  DESCRIPTIONS 

(C)  The  MSS  candidate  sensors  for  this  study  were  derived  from  the 
outputs  of  a pair  of  DIFAR  arrays.  Each  DIFAR  array  outputs  three 
channels  of  information.  The  first  channel  contains  the  output  of  an 
omnidirectional  hydrophone.  The  other  two  channels  contain  the  output 
of  collocated  horizontal  first-order  pressure  gradients.  One  gradient 
channel  has  its  main  response  axis  oriented  north/south,  and  the  other 
east/west.  The  output  on  each  of  these  channels  due  to  a plane  wave 
of  frequency  f and  amplitude  p which  is  incident  horizontally  to  the 
array  at  an  angle  6 east  of  north  can  be  expressed  as 

x (t)  = p siniot  Omnidirectional  Channel, 

o 

x^Ct)  = p cos6  sinut  North/South  (cosine)  Channel,  and 

x (t)  = p sinQ  sinwt  East/West  (sine)  Channel, 

s 

where  u=2irf.  These  gradient  channels  have  the  familiar  dipole 
directivity  pattern  (Fig.  1-3). 

(C)  It  will  be  assumed  that  the  response  function  of  each  channel  is 
rotationally  symmetric  about  its  main  response  axis.  Therefore,  if  the 
above  described  plane  wave  is  now  incident  at  an  angle  y above 
horizontal,  the  three  outputs  C3u  be  expressed  as 

x (t)  = p siuwt  , 
o 

x (t)  - p cos0  cosy  sinwt  , and 
c 

x (t)  ■ p sin0  cosy  sintot 
s 

It  should  be  noted  that  both  gradient  channels  null  out  any  signals 
arriving  from  directly  above  or  below  (y=±9Q0)  the  array. 

(U)  DIFAR  horizontal  beamforming  is  usually  accomplished  as  a weighted 
summation  of  the  three  channels: 
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(C) 


y (t,$)  “ [ax  (t)+cos<j)x  (t)+sin<{>x  (t)]/(l+a) 
o c 8 

= p sir.ut[ah:osY(cos4>  co&0-l-sin<p  sin0)]/(l+a) 


(1) 


The  present  formulation  assumes  that  both  directional  channels  are 
weighted  by  the  same  nondirectional  constant  which  has  been  equivalently 
represented  as  a single  omni  weight  factor  a.  When  a beam  is  formed 
towards  the  incoming  plane  wave  (4>=*Q ) , 

y(t,0)  = p sinu)t[ct+cosy]/(l+a) 


Note  that,  if  the  plane  wave  is  not  norizontaliy  incident  (y**0)»  the 
beamformer  response  will  be  less  than  chat  of  the  omnidirectional 
channel.  The  beam  patterns  of  the  above  described  DIFAR  beamformer  are 
termed  limacons.  The  shape  of  these  limacons  can  be  controlled  by  the 
omni  weight  factor  a.  Limacon  plane  wave  beam  patterns  for  typical 
omni  weight  factors  are  shown  in  Fig.  1-4.  Each  of  these  responses 
is  also  rotationally  symmetric  about  this  main  response  axis.  The 
limacon  formed  by  setting  a =>1 .0  is  referred  to  as  a cardioid. 

(C)  Some  of  the  MSS  candidate  sensors  currently  under  consideration 

consist  of  four  orthogonal  limacons.  Their  beamformer  can  be  expressed 
from  Eq.  (1)  as 


y (t ,o°)  » 

[axo(t)+xc(t}}/(l+a> 

North , 

y(t,9C°)  - 

[ax  (t)+x  (t)  l/d'Kx) 
o s 

Last , 

y(t,180°)  = 

[axo(cl-xc(t)j/a+a) 

South, 

y(t,270°)  - 

[axo(t)-xa{t)]/<Hn) 

West 

y 

CONFIDENTIAL 


1*0* 


a s 0.50  * * ■ 
a s 1.00  — 


FIGURE  1-4 

LIMACON  BEAM  PATTERNS 
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(C)  A simplified  block  diagram  of  this  beamformer,  through  Automatic  Line 
Integration  (ALI),  is  shown  as  Fig.  1-5.  When  the  omni  weight  factor 
is  set  to  1.0,  the  beamformer  output  consists  of  four  cardioids 
(Fig.  1-6).  This  configuration  is  termed  the  single  cardioids  sensor 
and  will  be  referred  to  throughout  this  report. 

(C)  The  time  domain  beamformer  discussed  earlier  is  not  an  efficient 

technique,  because  each  of  the  four  output  channels  must  be  Fourier 
transformed.  A more  efficient  technique  would  be  to  transform  each 
of  the  three  input  channels  before  beamforming.  Let  the  Fourier 
transform  of  each  DIFAR  channel  be  approximated  by 

N-l 

X(f)  “ h H x(t)e"iut/N 

t=0 

where  i“/^T.  Letting  Y(f)  be  an  observation  vector, 

-xo(f)- 

Y(f)  = Xc(f) 

X (f) 

L s J 

and  V ( tp ) be  a steering  vector, 


Much  of  the  following  discussion  follows  that  of  Ref.  2. 
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(C)  The  beamformer  output  becomes 

A(f,«)  “ V* (<f>)  Y(f) 

where  * denotes  the  complex  conjugate  transpose  operation.  For  sensors 
composed  of  our  orthogonal  limacons,  A( f , <J> ) becomes 


A(f ,0°)  = [aX  (f)+X  (f)] 
o c 

North 

A(f ,90°)  = [aX  (f)+X  (f)] 
o s 

East 

A(f ,180°)  = [aX  (f)-X  (f ) ] 

South 

o c 

A(f ,270°)  = [aX  (f )-X  (f ) ] 

West 

The  beamformer  angular  power  response  is  then  computed  as 

o(f,*)  = A(f  »$)A*(f  ,<j>) 

A block  diagram  of  this  Fourier  coefficient  beamformer  is  shown  as 
Fig.  1-7. 

(C)  The  ARL:UT  measurement  is  intended  to  serve  as  a tool  for 

Investigating  sensor  performance.  As  an  investigative  tool,  it  is 
desirable  to  save  the  DIFAR  data  at  a point  prior  to  where  it  is  formed 
into  specific  beams,  but  after  it  has  incurred  most  of  its  computational 
overhead  (i.e.,  FFl) , For  the  time  domain  beamformer,  this  would  require 
saving  the  time  series.  For  Fourier  coefficient  beamforming,  this  would 
require  saving  each  complex  Fourier  transform.  An  alternative  technique, 
which  meets  the  desired  criteria,  is  to  compute  and  save  cross-spectra 
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(C)  The  properties  of  DIFAR  cross-spectra  are  discussed  in  Ref.  2.  A few 
of  the  most  important  properties  are  summarized  below. 

(1)  Because  the  DIFAR  elements  are  spaced  much  less  than  a 
wavelength  apart,  the  acoustic  signals  on  all  three  channels  are  in  phase. 

(2)  Because  the  channels  are  all  in  phase,  each  element  q of 
Q is  real. 

(3)  Q is  symmetric  about  its  main  diagonal,  that  is, 

(4)  If  the  signal  and  noise  are  zero  mean  Gaussian,  then  0 contains 
all  of  the  information  contained  in  the  sequence  of  observations.  The 
DIFAR  beamformer  angular  response  now  becomes 


a(f,$)  « V*($)  Q(f ) V($) 


Carrying  out  the  indicated  matrix  operations  yields 


o(f,$)  * (C^+C2cos^+C^sin^+C^cos2^+C^sin2<>)/ (1+a)* 


(2) 


where 


C1  ■ “ ’ll  + 


C2  ■ 2oqu 

S • 2a<1n 

°4  • <q22-q}3)/2 

S ' q23 


It  should  be  noted  that  the  above  expression  has  the  form  of  a Fourier 
series  in  $.  This  implies  that  Q contains  information  about  an  isotropic 
terra,  as  well  as  two  directional  terms.  This  further  implies  that  this 
bearafortser  can  discriminate  between  two  signal  arrivals  in  the  same 
frequency  cell.  The  Fourier  series  also  suggests  that  an  inverse  FFT 
can  be  used  to  efficiently  forra  a large  number  of  equispaced  beams. 
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(C)  For  the  special  case  of  four  orthogonal  limacons,  the  beamformer 

outputs  become 

0(f,O°)  * cx  + c2  + c4  " (a2q11+2aq12'Kl22)/(l-»a)i  North  , 

o(f»90°)  * * (a2q11+2aq13+q33>/(l-ta)2  F.ast  , 

o (f ,180°)  “ Cx  - C2  + C4  - (a2qir2ctq12+q22)/(l4n)2  south,  and 

°(f»270°)  “ " c3  ~ c4  “ (c‘2q11-2aq13+q33)/(l+a)2  West 

A block  diagram  of  this  beamformer  is  shown  as  Fig.  1-8.  It  should  be 
noted  that  both  the  FFT  and  ALI  occur  prior  to  beam  formation.  This  is 
the  beamformer  configuration  used  for  this  study.  The  cross-spectra 
were  preserved  following  ALI.  Preservation  of  Q required  saving  six 
real  numbers  for  each  frequency  cell,  only  five  of  which  are  used  by  the 
four-1 imacon  beamformer.  The  complete  presentation  of  Q allows  beams  to 
be  formed  in  any  direction  via  Eq.  (2).  It  also  permits  application  of 
optimal  and/or  adaptive  techniques  such  as  those  described  in  Ref.  3 
and  applied  to  DIFAR  in  Ref.  2. 

(C)  The  various  1 imacon  beams  described  above  can  be  formed  from  a 

single  DIFAR  array.  Because  of  their  rotational  symmetry,  the  horizontal 
response  of  such  beams  is  also  their  vertical  response.  As  stated  earlier 
the  array  configuration  for  this  study  consisted  of  a pair  of  DIFAR  arrays 
These  arrays  were  arranged  one  above  the  other  with  a spacing  of  6.1  m 
(20  ft).  The  lower  DIFAR  array  (B)  was  used  as  described  above  to  form 
various  lloacons,  and  its  omnidirectional  channel  was  used  as  the 
reference  for  measuring  relative  performance  between  sensors.  Each 
channel  of  the  upper  DIFAR  array  (A)  was  subtracted  from  the  corresponding 
channel  of  the  lower  DIFAR  array  (B-A)  to  yield  what  is  termed  differenced 
DIFAR  channels.  This  can  be  expressed  as 
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(C)  xB  A(t)  * xV)  - xA(t)  Differenced  Omnidirectional  Channel, 
o o o 

XB  A(t)  » xB(t)  - xA(t)  Differenced  North/South  (cosine)  Channel,  and 
c c c 

XB  A(t)  ■ xB(t)  - xA(t)  Differenced  East/West  (sine)  Channel 
c s s 

The  purpose  of  this  arrangement  was  to  achieve  a horizontal  null  at  all 
frequencies.  The  horizontal  response  of  these  differenced  channels  is 
the  same  as  the  response  of  the  single  DIFAR  channels.  However,  the 
vertical  response  is  different  and  is  frequency  dependent. 

(C)  Given  a plane  wave  of  frequency  f incident  to  the  array  at  a 

horizontal  bearing  angle  0°  east  of  north,  and  a vertical  angle  Y°  above 
horizontal,  the  response  of  a single  DIFAR  array  llmacon  with  a horizontal 
steer  angle  of  $,  relative  to  the  response  of  the  omnidirectional  sensor, 
is 


b (0,Y,$,o)  • (a+cosY(eos$cos0+sin*sin$)]/(l-Ki)  . 

s 

The  angle  0 then  corresponds  to  the  horizontal  bearing  of  the  incident 
plane  wave  iu  degrees  east  of  north,  while  is  the  horizontal  steering 
angle  for  the  DIFAR  beam,  also  in  degrees  east  of  north.  The  angle  Y 
is  the  vertical  incidence  angle  of  the  plane  wave  measured  in  degrees 
above  horizontal.  Without  loss  of  generality,  ehoose  $-0*;  then 

b (9,Y»o)  * (a+cosi cos®)/(i+a) 
s 

The  response  of  the  differenced  omnidirectional  channel  is 

byyff.'t)  « sin^p^  siny)  , 

where  d is  the  hydrophene  separation  and  cq  is  the  sound  speed  in  the 
vicinity  of  the  array.  This  sensor  will  be  referred  to  throughout  this 
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(C)  report  as  the  vertical  dipole  sensor.  Vertical  plane  wave  beam 

patterns  at  approximate  frequencies  of  Interest  are  shown  in  Fig.  1-9. 
In  this  figure  and  Fig.  1-10,  0°  and  180°  correspond  to  horizontal 
arrival  angles,  while  90°  corresponds  to  a vertical  arrival  angle 
from  the  surface.  The  response  of  the  differenced  DIFAK  array  limacons 
can  bo  obtained,  by  the  product  theorem  (Ref.  4),  as  the  product  of  the 
single  limacon  and  vertical  dipole  responses, 

bD(f,9»Y,a)  * bs(0,y,a)  b^Cf.y) 

C J 

- (o+cosYcos0)  sin  - — sinY)/(l-K*) 

o 

The  only  differenced  limacons  included  in  this  study  were  differenced 
cardioids  (a*l).  Differenced  cardloid  vertical  plane  wave  patterns  at 
approximate  frequencies  of  interest  are  shown  in  Fig.  1-10. 

(C)  The  receiving  directivity  index  (Dl)  of  a beam  i3  defined  as 

(Ref,  4) 


DI  - 10  log 


/ 


fsM 

b2(0,Y>  da 


For  3-di©ensional  isotropic  noise,  the  single  limacon  directivity  index 
is  (Fig.  l-ll(a)) 

BI(a)  = 10  log  [{14n>y(a24~)] 

The  directivity  index  for  single  cardioids  (a*l)  is  4.77  dB.  The 
maximum  DI  occurs  at  a»l/3,  where  the  DI  is  6.02  dB.  Single  DtFAR  array 
limacons  formed  with  a*l/3  are  referred  to  in  this  study  as  maximum 
gain  limacons. 
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(C)  In  an  environment  where  all  of  the  noise  is  arriving  horizontally 

at  the  array,  it  is  appropriate  to  compute  DI  with  a 2-dimensional 
horizontally  isotropic  noise  field.  This  result  is  (Fig.  I-ll(a)) 

DI(cO  = 10  log^(14a)y^a2  + ijj  . 

The  corresponding  DI  for  single  cardioids  is  4.26  dB,  and  the  maximum 
DI  occurs  at  a=0.5,  where  the  DI  is  4.77  dB.  This  type  of  noise  field 
has  been  assumed  for  BEARTRAP  DIFAR/KANDE  processing  (Ref.  5). 

(C)  All  of  the  DIFAR  derived  sensors  considered  in  this  study  consist 

of  four  orthogonal  limacons.  The  minimum  horizontal  response  of  such 
sensors  will  occur  halfway  between  the  main  axis  of  adjacent  beams  (45°). 
The  response  at  selected  angles  relative  to  that  of  the  omnidirectional 
sensor  is  shown  in  Fig.  I-ll(b).  The  minimum  response  for  the  single 
cardioids  sensor  is  -1.38  dB,  and  is  -2.15  dB  for  the  maximum  gain  limacons 
sensor. 

(C)  The  directivity  index  of  the  vertical  dipole  sensor  in  a 

3-dimensional  isotropic  noise  field  is  (Ref.  6) 

Di(f)  = io 

L 0 O •* 

= 3.01  - 10  log[l-Jo(0.025f)) 

where  J (x)  is  the  Bessel  function  of  the  first  kind  and  zero  order, 
o 

(Recall  that  at  frequencies  below  approximately  125  Hz,  the  peak  array 
response  is  less  than  unity,  and  therefore  the  above  expression  is  not 
the  usual  DI.)  This  result  is  plotted  as  a function  of  frequency  in 
Fig.  I-ll(c).  This  sensor  was  primarily  intended  for  use  only  below 
2‘  Hz,  where  the  hydrophone  spacing  is  less  than  a wavelength.  At 
frequencies  less  than  approximately  125  Hz,  the  main  response  axis  is 
vertical  and  the  maximum  plane  wave  response  is  always  less  than  that 
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(C)  of  the  omnidirectional  sensor.  At  very  low  frequencies  (<10  Hz)  the 
hydrophone  spacing  is  only  a small  fraction  of  a wavelength;  thereto.?, 
the  vertical  dipole  response  to  plane  wave  inputs  should  be  very  small 
(<-18  dB//omni). 

(C)  The  directivity  index  of  the  differenced  limacons  is  (Ref.  6) 


DI(f,a)  = 10  log- 


4 ( 1-fa ) ‘ 


(2a2+l) 


/ sin2a\ 

-1  + 

[ 2*  'A 

\ 2a  ) 

3 

' 4a3  ' 

sin2a  + 


cos2a 
2a2  ' 


where  a=iTfd/c  — 0.0126f.  This  result  has  been  plotted  as  a function  of 
o 

frequency  (Fig.  I-ll(c))  and  omni  weight  factor  (Fig.  I-ll(d)).  The 
above  comments  on  the  frequency  dependence  of  the  vertical  dipole  also 
apply  to  the  differenced  limacons. 


(C)  The  actual  array  gains  of  the  single  DIFAR,  vertical  dipole,  and 

differenced  DIFAR  arrays  is  a function  not  only  of  the  directional 
properties  of  the  signal  and  ambient  sound  fields,  but  also  of  the 
interchannel  phase  and  amplitude  balance,  and  the  self-noise  of  the 
array.  A discussion  of  the  DIFAR  array  gain  degradation  due  to  these 
effects  is  found  in  Ref.  7. 


(S)  Following  signal  detection,  one  of  the  next  operations  to  he 

performed  for  each  DIFAR  sensor  is  bearing  estimation.  As  stated  earlier, 
the  current  MSS  approach  is  to  compose  each  DIFAR  sensor  of  four  orthogonal 
limacons.  Consider  now  a plane  wave  of  frequency  f incident  to  the 
single  DIFAR  array  at  a horizontal  angle  8 east  of  north  and  a vertical 
angle  y above  horizontal.  In  the  absence  of  noise,  the  four  beamformer 
outputs  will  be 
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(S) 


o(f,0#) 
a(f ,90°) 
0 (f ,180°) 

o(f,270°) 


(a+cosYC088  )^  / (14a)^ 
2 2 

(a+cosYsinO)  /(1+a) 

2 2 

(o-co8YCO80)  / (1+a) 
(a-cosYsinO)^/  (].+a)^ 


North  , 
East  , 
South,  and 

West 


The  current  bearing  estimator  is  (Ref.  1) 


B » tan 


■ tan 


-lfo(f.90°)  - a (f ,270° ) 
L o (f ,0°)  - o(f,180°) 

-if  2cosYsin8  "I 
L 2COSYCO30  J 


] 


« 8 


Real  signals,  however,  are  seldom  encountered  without  any  noise.  To 
minimize  the  effect  of  noise  on  the  estimated  bearing,  an  estimate  of 
the  background  noise  level  on  each  beam  is  subtracted  from  the  beamformer 
outputs  prior  to  bearing  estimation 

og(f,$)  - os+N(f,^)  " 

where  the  bar  indicates  an  estimated  quantity.  This  is  sometimes  termed 
acoustic  debiasing.  The  same  bearing  estimator  can  be  shown  to  apply  to 
the  differenced  limacon  sensors. 

(C)  As  stated  earlier,  the  DIFAR  cross-spectral  beamformer  can 

discriminate  between  two  signal  arrivals  in  the  same  frequency  cell.  A 
maximum  entropy  bearing  estimation  technique  for  DIFAR  data  has  been 
developed  (Ref.  5)  which  estimates  two  bearings  directly  from  the 
cross-spectral  matrix.  This  technique  could  be  Investigated  for  applica- 
tion to  KSS.  Another  technique  has  been  developed  by  Pryor  (Ref.  8), 
where  « generalized  limacon  is  used  to  automatically  null  out  the  effects 
of  broadband  directional  noise.  When  the  omnidirectional  channel  is 


26 


CONFIDENTIAL 


(C)  replaced  with  this  generalized  limacon,  the  bearings  re  automatically 
debiased  and  detection  performance  improved.  Thiv  technique  may  also 
merit  further  investigation. 

(U)  The  principal  characteristics  of  each  of  the  sensors  included  it. 

this  study  are  given  in  Table  1-1. 

(C)  In  the  process  of  documenting  this  study,  it  became  apparent  that 

a sensor  with  a response  similar  to  that  of  the  vertical  dipole  could 
be  formed  from  a single  DIFAR  array.  Consider  a beam  formed  from  the 
power  on  each  of  the  three  DIFAR  channels  by 

X*(f)  = X*(f)  - X*(f)  - X*(f) 

Substituting  the  response  of  each  of  the  DIFAR  channels  yields 

X*(f)  - 1 - cos^Q  cos^y  - sin^B  cos^y 

, 2 
■ 1 - cos  y 

, 2 
» sin  Y 

Thus  the  response  of  this  sensor  is 

by(Y)  * sin j y I , 

where  y is  the  vertical  angle  above  horizontal. 

(C)  This  sensor  has  several  advantages  over  the  vertical  dipole:  Because 

only  power  terms  are  used,  it  is  insensitive  to  phase  imbalance  between 
DIFAR  channels;  it  has  the  same  response  at  all  frequencies;  and  it  does 
not  rely  on  the  coherence  of  the  acoustic  field  across  a 6.1  a (20  ft) 
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Differenced  Cardioids  4-orthogonal  differenced  (B-A)  Yes  10  to  250  14.34  at  50  Hz 

(DC)  cardioids  (a“=1.00)  providing  a 7.43  at  150  Hz 

horizontal  null  at  all  frequencies  8.21  at  300  Hz 
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(C)  aperture.  The  DI  of  this  sensor  is  the  same  as  that  of  a gradient 
channel  (4.77  dB).  Because  of  time  constraints,  this  sensor  was  not 
included  in  this  study,  but  does  appear  to  merit  further  investigation. 
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III.  DATA  COLLECTION  AND  PREPROCESSING 

(C)  Data  collection  was  accomplished  via  three  specially  modified 

ACODAC  (Acoustic  Data  Capsule)  systems.  Each  ACODAC  system  consisted 
of  an  instrumentation  pressure  vessel  (IPV)  containing  a 7-channel 
tape  recorder  and  associated  signal  conditioning,  timing,  and  control 
circuits  (Ref.  9).  Prior  to  the  FVT  exercise,  extensive  modifications 
were  made  to  these  systems  (Ref.  10).  The  resultant  array  configura- 
tion is  shown  in  Fig.  1-12.  DIFAR  arrays  were  located  18.3  m (60  ft) 
and  24.4  m (80  ft)  above  the  ocean  floor.  The  electronics  were 
contained  within  the  sensor  housings,  junction  boxes,  and  IPV  (Fig.  1-13). 
The  directional  channels  of  each  DIFAR  array  were  independently  compass 
corrected  to  magnetic  north.  Each  of  the  six  alog  data  channels  was 
summed  with  an  amplitude  modulated  reference  tone.  The  seventh  channel 
contained  IRIG  C timecode  which  encoded  the  day,  hour,  minute,  system  ID, 
and  AGC  settings.  The  data  collection  characteristics  are  summarized 
in  Table  1-2. 

(C)  The  two  pressure  gradient  channels  of  each  DIFAR  array  were  each 

obtained  from  the  difference  of  two  geophones.  Since  these  are  velocity 
sensors,  they  will  be  even  more  sensitive  to  flow  noise  than  hydrophones. 

(C)  Because  of  the  very  low  ACODAC  tape  speeds  (normally  15/160  ips), 

significant  crosschannel  alignment  errors  are  incurred  during  recording, 
duplication,  and  reproduction  of  these  tapes.  Since  the  FVT  ACODAC  data 
were  primarily  intended  for  DIFAR  processing,  special  procedures  were 
required  to  improve  crosschannel  alignment.  The  calibration  accuracy 
required  to  achieve  acceptable  DIFAR  performance  (Ref.  10)  is  summarized 
in  Table  1-3.  Achieving  this  accuracy  from  the  ACODAC  recordings  required 
a special  playback  facility.  This  playback  facility  used  the  reference 
signal  on  each  data  channel  to  align  the  three  channels  of  each  DIFAR 
set.  This  procedure  has  frequently  been  termed  D"*  (dewow,  deskew,  and 
deflutter).  The  aligned  outputs  of  this  facility  were  re-recorded  in 
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(C)  TABLE  1-2 

DATA  COLLECTION  SYSTEM  CHAFACTERISTICS  (U) 


DATA  BAND: 

10  to  600  Hz 

DYNAMIC  RANGE: 

60  dB  (Recorder  27  dB,  AGC  30  dB) 

RECORDER  SPEED: 

15/80  ips 

RECORDING  MODE: 

Analog  Direct 

RECORDING  DURATION: 

6.5  days  (maximum) 

AGC: 

5 to  35  dB  in  10  dB  step  with  3 dB  hysteresis  and 
gain  of  direct ionals  slaved  to  omnis 

ARRAY  GEOMETRY: 

Lower  DIFAR  (B)  18.3  m (60  ft)  above  bottom 

Upper  DIFAR  (A)  24.4  m (80  ft)  above  bottom 

Pressure  Vessel  (IPV)  36.6  m (120  f ) above  bottom 

NUMBER  OF  CHANNELS: 

7 (1/2  in.  tape) 

CHANNEL  ASSIGNMENTS: 

1 B cos  (North/South) 

2 B-A  sin  (East/West) 

3 B omni 

4 B-A  omni 

5 B sin 

6 B-A  cos 

7 IRIG  C time  code  including  AGC  settings 
(only  AGC  settings  for  channels  3 and  4 
are  meaningful) 

REFERENCE  SIGNAL: 

1-50  Hz  sine  wave  with  125  Hz  amplitude  modulation 
recorded  in  parallel  on  all  data  channels  (1 
through  6) 
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(C)  TABLE  1-3 

SENSOR  CALIBRATION  ACCURACY  REQUIREMENTS  (U) 


SINGLE  DIFAR  ARRAY  - 

sin  to  cos  Tracking: 
sin  to  ooni  Tracking: 
Compass  Trim: 

Frequency  Range: 

DIFFERENCED  DIFAR  ARRAY  - 

omni  to  omni  Tracking: 
cos  to  cos  Tracking: 
sin  to  sin  Tracking: 
Tilt: 


15°  phase  ±1  dB 
20°  phase  ±3  dB 
5°  p-p  phase  ±1  dB 
20  to  600  Hz 


15°  phase  +0.75  dB 
7.5°  phase  +0.5  dB 
7.5°  phase  +0.5  dB 
~3°  maximum 
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(C)  FM  mode  at  a much  higher  tape  speed.  Recordings  of  this  type  are 
relatively  immune  from  phase  and  amplitude  errors.  As  many  as  20  FM 
tapes  were  produced  from  each  ACODAC  tape. 

(C)  After  a single  set  of  aligned  tapes  was  produced,  the  tapes  were 

duplicated  and  shipped  to  the  various  data  processors,  including 
ARL:UT.  The  duplication  of  FM  tapes  is  usually  a simple  procedure. 

The  duplication  station  is  usually  configured  to  perform  a unity-gain 
DIRECT-TO-DIRECT  mode  copy.  Since  the  FM  signal  is  not  decoded,  no 
additional  errors  are  incurred.  Unfortunately,  the  duplication  station 
used  for  these  data  did  not  have  the  required  DIRECT  mode  capacity; 
therefore,  FM-to-FM  copies  were  made.  The  resultant  copies  were  of 
poor  quality.  Small  portions  of  the  data  were  completely  destroyed  due 
to  overdeviation  of  the  FM  carrier.  Large  amplitude  imbalances  (6  to 
8 dB)  were  incurred,  even  between  two  copies  of  the  same  tape.  Conse- 
quently, the  duplication  procedure  may  have  significantly  increased  the 
amplitude  error  of  these  results,  which  would,  in  turn,  have  increased 
bearing  error. 
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IV.  A/D  CONVERSION  AND  DATA  REDUCTION 


(U)  The  first  major  component  of  the  ARL:UT  measurement  system  (Fig.  1-2) 

is  termed  A/D  conversion  and  data  reduction.  This  component  accepts 
single  and/or  multichannel  inputs  in  analog  format  and  converts  them  to 
time-averaged  calibrated  autospectra  and  cross-spectra-  These  spectra 
serve  as  the  primary  input  to  both  of  the  remaining  measurement  system 
components.  For  this  application,  the  required  steps  are  shown  in 
Fig.  1-14. 

(S)  Selected  portions  of  the  FVT  ACODAC  data  set  were  processed  with 

two  frequency  resolutions.  The  low  resolution  data  were  intended  to 
emulate  MSS  standard  resolution  (0.2  Hz)  data.  The  high  resolution 
data  were  intended  to  emulate  MSS  vernier  resolution  (0.0125  Hz)  data. 

(C)  The  selected  data  were  first  converted  from  an  analog  to  a digital 

representat  ".on  and  stored  as  an  intermediate  data  product  on  magnetic 
tape.  The  available  timecode  md  calii'ration  information  were  simultan- 
eously extracted.  The  analog  data  were  reproduced  on  an  Ampex  instru- 
mentation tape  recorder  with  intermediate  band  electronics.  Three 
channels,  a DIFAR  set,  were  digitized  on  each  pass  through  the  data 
set.  The  analog  hardware  configuration  was  dependent  on  whether  standard 
frequency  resolution  (Fig.  1-15)  or  vernier  frequency  resolution 
(Fig.  1-16)  was  desired.  The  analog  signals  were  bandpass  filtered  to 
prevent  aliasing  and  input  to  programmable  amplifiers  under  joint 
computer/operator  control.  The  amplifier  gain  was  automatically  varied  at 
1 rain  intervals  to  obtain  maximum  utilization  of  the  dynamic  range  of  a 
12-bit  A/D  converter.  The  gains  of  the  gradient  channels  were  slaved  to 
the  gain  of  the  omnidirectional  channel,  except  for  constant  offsets  which 
were  frequently  required  to  bring  the  channels  into  approximate  balance.  The 
A/D  converter  was  preceded  by  three  simultaneous  sample  and  hold  (SS-t-H) 
circuits,  one  per  channel.  Vernier  resolution  data  were  also  frequency 
translated  prior  to  the  SS+H  circuits.  The  frequency  translator  consisted 
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FIGURE  1-14 

A/D  CONVERSION  ANO  DATA  REDUCTION  OVERVIEW 
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(C)  of  three  custom  made  double-balance  mixers  and  low  pass  filters 

(Ref.  11).  The  gain  of  these  low  pass  filters  was  adjusted  so  that 
the  linear  input  range  of  the  mixer  corresponded  to  the  range  of  the 
A/D  converter.  The  reference  for  these  mixers  was  derived  from  the 
timebase  reference  signal  recorded  on  each  analog  tape.  The  SS+H 
circuits  and  A/D  converter  were  also  triggered  at  a rate  derived  from 
this  reference. 

(U)  The  digitized  data  were  then  input  to  a Control  Data  3200  computer 

system  (Fig.  1-17).  This  computer  was  used  to  select  amplifier  settings 
and  to  ensure  that  the  data  did  not  exceed  the  ±10  V range  of  the  A/D 
converter.  The  operator  had  the  option  of  permitting  such  "clipping," 
and  usually  did  when  high  level  transient  signals  were  present.  Such 
"clipped"  data  were  automatically  tagged  and  then  deleted  during  data 
reduction.  The  digital  data  were  stored  on  magneLic  tape  in  a format 
designed  to  facilitate  later  processing,  along  with  time,  gain,  and 
quality  indicators.  The  computer  performed  error  detection  which 
included  timecode  monitoring  and  verification,  signal  amplitude  monitor- 
ing, data  flow  monitoring,  and  overall  synchronization.  If  synchro- 
nization was  lost,  the  process  was  terminated,  and  the  operator  was 
informed  of  its  cause.  After  the  appropriate  corrective  action  was 
performed,  the  process  was  reinitiated  prior  to  the  point  where 
synchronization  had  been  lost.  The  computer  also  generated  a permanent 
log  of  each  interruption  and  the  operator's  reaction. 

(U)  The  data  were  digitized  in  a sequence  of  four  A/D  runs,  one  for  each 

of  the  three  vernier  resolution  frequency  bands  and  one  for  the  standard 
resolution  frequency  band.  The  playback  tape  recorder  was  first  equalized 
so  chat  the  header  calibration  sequences  had  a flat  response  across  the 
frequency  band  of  interest.  The  three  analog  data  channels,  from  tape 
recorder  output  to  SS+H  circuit  output,  were  then  balanced  to  within 
0.1  dB  and  ±3°.  The  desired  data  from  each  site  were  then  digitized. 

The  external  calibration  sequence  (header  cals)  for  these  data  were 
digitized  immediately  before  and  after  the  data  digitization.  This 
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(U)  redundance  was  used  to  verify  the  stability  and  repeatability  of  the 

playback  system.  Amplitude  differences  were  typically  less  than  0.1  dB. 

Interchannel  phase  differences  also  remained  nearly  constant  (Appendix 

3 

I-A) . The  amplitude  transfer  functions  (T)  of  the  entire  record-D  - 
duplicate-playback-A/D  system  were  obtained,  from  the  measured  levels  of 
the  calibration  signals  (M)  and  their  absolute  levels  (A) 


T (f ) 


A (f) 
M(f) 


> 


where  f is  frequency  in  hertz.  The  composite  transfer  function  (C) 
used  to  calibrate  the  data  to  sound  pressure  level  (SPL)  in  the  water 
was 


C(f) 


lii) 

R(f) 


* 


where  R(c)  is  the  amplitude  transfer  function  from  the  water  to  the 
tape  recorder,  exclusive  of  AGC.  A(f)  was  furnished  by  Woods  Hole 
Oceanographic  Institution,  and  R(f)  is  tabulated  in  Ref.  10. 

(C)  The  tabulated  values  of  R(f)  apply  only  to  single  DIFAR  arrays 

and  neglect  the  effect  of  the  differencing  network  (Fig.  1-13).  The 
calibration  of  the  differenced  DIFAR  data  (B-A)  was  accomplished  by 
using  the  linear  summation  of  the  corresponding  values  of  R(f). 

The  composite  calibration  factors  were  therefore  computed  by 


cB(f) 


yo 

Rb(0 


A(f) 

tM8(f)  RB(f)]  . and 
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(C)  „ _ TB-A(f) 

CB-A(f)  * [RA(f)+RB(£)] 

{Va(£)  ‘8A<£>+V£)'} 

As  was  discussed  in  section  II,  at  frequencies  of  less  than  125  Hz  the 
response  of  the  differenced  DIFAR  array  channels  will  always  be  less 


apart. 


(U)  Transfer  functions  for  the  standard  frequency  resolution  data  were 

computed  by  log-log  interpolation  on  T and  R.  The  number  of  calibration 
signals  within  each  vernier  frequency  band  was  inadequate  to  define  its 
response.  Therefore,  immediately  prior  to  and  following  the  digitiza- 
tion, the  tape  recorder  was  replaced  by  a white  noise,  generator,  and 
the  ARL:UT  A/D  system  response  was  measured.  These  curves  were  then 
used  to  interpolate  between  calibration  signals. 

(C)  The  final  data  reduction  step  accepts  as  input  the  digital  time-series 
data  and  composite  transfer  functions,  and  then  computes  the  outputs 
time-averaged  calibrated  cross-spectra.  All  of  the  computations  of  the 
cross-spectral  matrix  beaoformer  (Fig.  1-8)  through  ALI  are  performed. 

In  addition,  computations  are  performed  to  edit  and  normalize  the  input 
data  and  to  calibrate  the  results.  The  FFT  length  (8192  samples)  was 
chosen  in  conjunction  with  the  sample  rates  to  yield  the  desired 
frequency  range  and  resolution.  The  spectral  window  (Hanning)  and 
overlap  (50Z)  were  chosen  for  their  variance  reduction  properties 
(Refs.  12  and  13)  and  low  spectral  bias  (Ref.  14).  The  ALI  type  was 
rectangular  integration,  where  all  of  the  spectra  computed  during  the 
ALI  Interval  arc  averaged  across  each  frequency  cell.  The  spectra  output 
for  each  ALI  interval  arc  characterized  by  the  number  of  equivalent 
degrees  of  freedom.  These  parameters  are  summarized  in  Table  1-4. 
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(S) 


TABLE  1-4 


A/D  CONVERSION  AND  DATA  REDUCTION  PARAMETERS  (U) 


PARAMETER 

STANDARD  RESOLUTION 

VERNIER  RESOLUTION 

Sample  Rate 

1250  Hz 

100  Hz 

Frequency  Range 

10  to  600  Hz 

46  to  84  Hz  (low  band) 

145  to  183  Hz  (mid  band) 
300  to  338  Hz  (high  band) 

FFT  Length 

6.55  sec 

81.92  sec 

Spectral  Window 

Hanning 

Hanning 

Frequency  Spacing 

0.1526  Hz 

0.0122  Hz 

3 dB  Bandwidth 

0.2197  Hz 

0.0176  Hz 

Equivalent  Noise  Bandwidth 

0.2289  Hz 

0.0183  Hz 

FFT  Overlap 

50% 

50% 

ALI  Interval 

1 min 

5 min 

AH  Type 

Rectangular  Integration 

Rectangular  Integration 

FFT/ ALI  (maximum) 

18 

7 

Time  Bandwidth  Product  (maximum) 

13.67 

5.76 

Number  of  Equivalent 

32.3 

13.4 

Degrees  of  Freedom  (maximum) 


UNCLASSIFIED 


The  number  of  equivalent  degrees  of  freedom  of  each  spectra 
is  a measure  of  statistical  stability.  It  has  been  shown  (Ref.  15) 
that  overlapped  FFT  processing  of  windowed  Gaussian  data  results  in 
spectral  estimates  whose  distribution  is  very  closely  approximated  by 
a chi-square  distribution  with  the  appropriate  number  of  equivalent 
degrees  of  freedom.  Once  this  distribution  is  known,  it  can  be  used 
to  set  detection  thresholds  and  to  test  data  for  inhomogeneity.  The 
number  of  equivalent  degrees  of  freedom  for  n transforms  with  502 
overlap  and  Hanning  is  approximately  (Ref,  15) 


The  value  of  n may  vary  between  spectra  because  of  the  editing 
criteria  applied  to  the  input. 
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V.  AMBIENT  SOUND  FIELD  MEASUREMENT  TECHNIQUE 

(U)  The  second  major  component  of  the  ARL:UT  measurement  system 

(Fig.  1-2)  is  that  which  accepts  spectra  as  input,  and  computes  and 
displays  measurements  of  the  ambient  sound  field  (ASF).  The  goal  is 
to  obtain  measurements  as  a function  of  frequency  and  time  of  the  ASF 
levels  observed  by  the  omnidirectional  sensor  and  of  the  noise  gains 
of  the  directional  sensors  relative  to  the  omnidirectional  sensor. 
Clutter  measurements  are  also  a property  of  the  ASF;  however,  since 
such  measurements  require  a cw  processor,  discussion  of  them  will  be 
deferred  to  the  section  on  cw  processing.  Noise  gain  is  herein  defined 
as 


NG  ■ 10  loe  ASF  level  observed  via  directional  sensor 

ASF  level  observed  via  omnidirectional  sensor  ’ 

where  both  levels  are  determined  over  the  same  time  interval  and 
frequency  band.  If  the  ASF  was  isotropic,  noise  gain  would  be  equal 
to  the  negative  of  the  directivity  indices  computed  in  section  II: 


NG 


isotropic 


-DI 


(U)  The  ASF  measurements  were  formed  over  fixed  frequency  bands.  For 

the  vernier  resolution  data,  each  band  was  1 Hz  wide  and  centered  at 
each  integer  frequency  within  the  band.  For  the  standard  resolution 
data  each  band  was  1/10-octave  wide  (W»Q,1).  The  center  frequency 
(f  ),  lower  limit  (f  ),  and  upper  limit  (f  ) for  each  band  are  obtained 

C Ju  U 

from 
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(U) 


K - 2 


W/2 


fc(i)  - 10 


0.3W1 


fL(i)  = fc(D/K 


f (i)  - f (i)  x k 
u c 


where  i is  the  standard  band  number,  and  ranged  from  43  (fc  ■ 19.5  Hz)  to 

92  (f  ■=  575  Hz), 
c 


(C)  Techniques  for  extracting  a representative  cross-spectral  matrix 

from  a frequency  band  of  such  matrices  are  not  well  established.  If  such 
techniques  were  available,  then  the  resulting  matrix  could  be  used  to 
characterize  horizontal  directivity  and  optimum  beamforming  (Ref.  2). 

For  this  study,  the  elements  of  each  cross-spectral  matrix  were  summed 
across  the  desired  frequency  band;  the  resulting  sums  were  then  normalized 
to  a 1 Hz  bandwidth.  The  resulting  matrices  from  both  the  single  and 
differenced  DIFAR  arrays  were  then  merged  (Fig.  1-18)  into  a single  ASF 
data  base.  This  data  base  was  then  acted  upon  with  the  cross-spectral 
limacon  beamformer  (Fig.  1-8)  to  _ ield  beam  1_ -’els  and  noise  gains  as  a 
function  of  frequency,  time,  site,  censor,  direct. on,  and  omni  weight 
factor.  A description  of  the  three  i 'splay  foi  .a^ts  can  be  found  in 
Volumes  II  and  III. 


(U)  The  expected  distribution  about  the  mean  value  for  the  ASF 

measurements  through  l/10-i  ■'  ave  bands  for  stationary  Gaussian  input  data 
is  shown  i?i  Fig.  1-19.  The.-e  results  arc  derived  in  Appendix  I-B  and 
arc  based  on  Ref.  16.  Whenever  the  measured  distribution  exceeds  those 
of  Fig.  1-19,  this  implies  chat  the  input  data  were  nonstationary  or 
non-Gausslan  (cw  signals).  The  vertical  line  on  Fig.  1-19  denotes  the 
expected  spread  of  the  ernier  resolution  measurements. 
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FIGURE  1-18 

ASF  PROCESSOR  OVERVIEW 
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FIGURE  1-19 

EXPECTED  DISTRIBUTION  OF  AMBIENT  SOUND  FIELD 
MEASUREMENTS  THROUGH  1/10-OCTAVE  BANDS 


AS-79-41 
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(U)  An  additional  test  was  applied  to  these  data  to  ensure  that 

the  measurements  were  not  contaminated  by  the  electronic  noise  floor  of 

3 

the  record-D  -duplicate-playback-A/D  system.  The  level  of  the  electronic 
noise  floor  (E)  was  approximated  from  the  measured  levels  during  the  first 
minute  of  an  internal  calibration  signal  (6  h cal).  During  this  minute  an 
open  circuit  was  input  to  the  ACODAC  data  amplifiers.  This  electronic 
noise  is  always  present  and  biases  the  resultant  measurements  (R)  above  the 
actual  levels  (A), 


R = A + E 

If  biases  of  less  than  1 dB  only  are  to  be  tolerated, 

10  log(R/A)  < 1 

then  the  acceptable  measurement-to-noise  floor  ratio  is 

10  log(R/E)  < -10  logd-lO-0*1)  ■ 6.87  dB 

For  these  data,  any  measurements  within  6 dB  of  the  measured  electronic 
noise  floor  levels  were  deleted.  This  resulted  in  the  elimination  of  all 
of  the  data  below  30  Hz. 

(U)  The  ASF  processing  parameters  are  summarized  in  Table  1-3. 
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AMBIENT  SOUND 

Parameter 

ALI  Interval 

Bandwidth 

Band  Formation 
Technique 

Noise  Floor 

Contamination  Criteria 


TABLE  1-5 

FIELD  PROCESSOR  PARAMETERS 

Standard  Resolution 
1 min 

0.1  octave 

Rectangular 

Integration 

6 dB 
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(U) 

Vernier  Resolution 

5 min 
1 Hz 

Rectangular 

Integration 

6 dB 
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VI.  cw  MEASUREMENT  TECHNIQUE 

(U)  The  third  major  component  of  the  ARL:UT  measurement  system 

(Fig.  1-2)  accepts  spectra  as  input,  and  computes  and  displays 
measurements  based  on  received  cw  signals.  All  cw  signal  processors 
include  the  functions  of  signal  identification  and  parameter  estimation. 
The  selection  of  an  identification  technique  is  based  on  the  amount  of 
information  known  about  the  signal.  The  parameter  estimation  technique 
is  selected  to  provide  the  desired  measurements.  Different  processing 
techniques,  when  applied  to  the  same  data,  can  yield  significantly 
different  results. 

(C)  In  an  earlier  ARL:UT  FVT  study  (Ref.  17),  the  goal  of  the  cw 

processor  was  to  obtain  cw  propagation  loss  measurements.  It  was  assumed 
that  the  frequency  and  nominal  Doppler  shift  were  known  a priori  for  the 
single  projector  line  of  interest.  The  identification  technique  was 
a peak  search  over  a restricted  frequency  range,  and  explicit  line 
tracking  from  one  ALI  interval  to  the  next  was  not  required.  The 
propagation  loss  parameter  was  computed  from  the  cw  SPL,  which  was 
estimated  by  summing  the  received  SPL  over  the  band  of  frequency  cells 
containing  the  signal  and  subtracting  an  estimate  of  the  ASF  SPL  in  the 
same  band.  Measurements  obtained  with  this  technique  are  relatively 
insensitive  to  Doppler  smearing  and  to  the  distribution  of  signal  SPL 
between  and  within  each  cell.  The  error  of  such  measurements  is  further 
decreased  relative  to  these  based  on  single'  cells,  because  the  ASF  was 
summed  across  multiple  cells,  thus  decreasing  the  expected  error  of  a 
mean  estimate.  Finally,  the  measurements  based  on  total  received  signal 
SPL  correspond  to  those  predicted  by  most  propagation  models. 

(C)  In  this  study,  the  goal  was  to  obtain  performance  measurements  for 

candidate  MSS  sensors  coupled  to  a surveillance  processor.  This  processor 
had  no  a priori  information  on  tire  number  of  sources,  their  frequencies, 
or  their  locations.  It  was  required  to  detect  and  track  all  of  the  lines 
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within  the  processor  frequency  band.  The  desired  performance  measurements 
included  detection  and  bearing  estimation  against  simulated  target  lines, 
array  gain  of  candidate  sensors,  and  clutter  related  (processor  loading) 
statistics.  The  processor  was  also  required  to  output  sufficient 
ancillary  data  to  permit  an  interpretation  of  the  performance  measurements. 

These  cw  measurements  were  provided  by  a sequence  of  three  processors 
(Fig.  1-20).  The  single  cell  processor  determined  which  of  the  frequency 
cells  contained  signals,  and  estimated  their  parameters.  The  line 
related  processor  formed  the  detected  cells  into  lines,  estimated  line 
parameters,  and  tracked  the  lines  through  time.  The  editing  and  display 
processors  were  used  to  identify  projector  lines  and  exercise  artifacts, 
and  to  extract  and  display  the  final  cw  signal  and  clutter  measurements. 

Each  of  these  processors  was  implemented  as  a separate  software  program 
and  will  now  be  described  in  detail. 

Two  versions  of  the  single  cell  processor  were  used.  The  first 
version  provided  only  for  a single  beam  (omnidirectional  or  vertical  dipole) 
sensor  without  bearing  estimation,  while  the  second  was  a 4-beam  (single 
or  differenced  DIFAR)  sensor  version  with  bearing  (Fig.  1-21).  Only 
the  multibeam  version  will  be  discussed  since  it  contains  all  the 
features  of  the  single  beam  version. 

The  first  function  of  the  single  cell  processor,  beamforming,  was 
accomplished  with  the  cross-spectral  beamformer  discussed  in  section  II. 

The  SPL  values  for  each  beam  were  computed  in  parallel  (Fig.  1-8)  and 
presented  to  the  autodetect  routine  in  order  of  decreasing  frequency. 

The  mean  background  ASF  SPL  on  each  beam  was  estimated  by  a 
clipper/replacement  technique  (Fig.  1-22).  This  technique  synthesizes 
a two-pass  estimator.  On  the  first  pass,  a "clipped"  spectrum  (y)  is 
created  where  the  narrow  high  level  lines  have  been  replaced  with  lower 
levels.  On  the  second  pass,  the  mean  background  ASF  SPL  is  estimated  from 


54 

CONFIDENTIAL 


STANDARD  FREQUENCY  RESOLUTION 


CONFIDENTIAL 


CONFIDENTIAL 


FIGURE  1-22 

BACKGROUND  NOISE  ESTIMATION  AND  SIGNAL  DETECTION  ALGORITHM  (U) 


CONFIDENTIAL 


(C)  the  clipped  spectrum.  This  technique  was  intended  to  reduce  the  ASF 
estimator  bias  caused  by  high  level  lines. 

(C)  During  the  first  pass  of  the  estimator,  the  unclipped  spectrum  (x) 

was  first  averaged  over  a 32-cell  window  centered  about  each  test  cell, 

i+15 


If  the  ratio  of  the  test  cell  SPL  to  the  average  unclipped  SPL  exceeded 
a threshold,  termed  the  clipper  coefficient, 


9 


then  the  test  cell  SPL  was  replaced  by  the  product  of  the  average 
unclipped  SPL  and  the  replacement  coefficient. 


otherwise 


(C)  During  the  second  pass  of  the  estimator,  the  background  ASF  estimate 

was  obtained  by  averaging  the  clipped  spectrum  over  a window  of  32  cells 
for  standard  resolution  and  64  cells  for  vernier  resolution.  Each 
window  was  centered  about  a 3-cell  gap  which  contained  the  detection 
test  cell. 

(C)  Selection  of  clipper  and  replacement  coefficients  requires  a tradeoff 

in  the  accuracy  of  the  estimator  against  high  level  and  low  level  lines. 
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(C)  By  reverting  once  more  to  the  chi-square  approximation  for  FFT 

processing  of  stationary  Gaussian  input,  the  clipper  and  replacement 
coefficients  can  be  computed  as  functions  of  ASF  percentile  level  and 
the  number  of  equivalent  degrees  of  freedom  (Fig.  1-23).  With  a 
clipper  probability  of  0.90,  10%  of  the  test  cells  would  be  clipped  in 
the  absence  of  any  signals.  This  clipping  of  nonsignal  cells  would 
tend  to  bias  the  ASF  estimate  low  (Appendix  I-C).  Any  decrease  in  the 
clipper  probability  would  increase  this  negative  bias  (Fig.  1-24). 

The  effect  of  a low  ASF  estimate,  in  the  absence  of  signals,  is  to 
increase  the  false  alarm  rate  since  the  detection  decision  is  based 
upon  the  estimated  ASF  levels.  Choosing  a higher  clipper  probability 
would  decrease  the  negative  bias  in  the  absence  of  signals,  but  would 
also  tend  to  increase  the  positive  bias  of  the  estimator  when  signals 
are  present. 

(C)  This  estimator  performs  well  against  very  narrow  signals,  but 

degrades  quickly  if  the  signal  occupies  more  than  a few  cells.  This 
degradation  results  in  high  ASF  estimates  in  the  vicinity  of  broad 
and/or  very  high  level  signals.  The  result  of  this  positive  bias  is 
to  prevent  detection  of  weak  signals  near  strong  signals  and  to  cause 
errors  in  the  estimated  signal  parameters  (Ref.  18).  These  biases 
will  be  larger  for  a square-law  detector,  such  as  was  used  for  this 
study,  than  for  a linear  detector.  For  this  study  the  clipper 
probability  was  0.90,  the  replacement  probability  was  0.50,  and  the 
detection  coefficient  was  adjusted  to  compensate  for  the  expected 
estimator  bias  (Fig.  1-24). 

(C)  The  difficulty  of  automatically  obtaining  an  unbiased  background 

ASF  estimate  is  well  known.  It  has  been  shown  (Ref.  18)  that  a 
median  based  estimate  will  have  a ouch  smaller  bias;  however,  the  normal 
computation  teciwiquc  of  ordering  the  cells  within  the  noise  band  is 
too  costly  for  operational  realtime  systems.  An  alternative  median 
estimation  technique,  which  is  less  costly  than  even  the  clipper/ 
replacement  technique,  has  recently  been  suggested.  For  a reasonably 
flat  spectrum,  the  suggested  procedure,  following  initialization,  is  to 
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FIGURE  1-23 

CLIPPER  AND  REPLACEMENT  COEFFICIENTS  VS  NUMBER  OF 
DEGREES  OF  FREEDOM  RT  SELECTED  PROBABILITIES  (U) 
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FIGURE  1-24 

BIAS  IN  ESTIMATED  NOISE  LEVEL  DUE  TO  CLIPPING 
OF  SPECTRUM  AT  VARIOUS  PROBABILITY  LEVELS 
WITH  MEDIAN  VALUES  AS  REPLACEMENTS 
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(C)  sequentially  compare  each  cell  of  the  undipped  spectrum  (x)  with 
the  current  value  of  the  estimator  (X).  Depending  upon  the  result  of 
this  comparison,  the  estimator  is  updated  as  follows: 

if  > X^_^,  then  ■ X^_^  + a ; 
if  t*ien  * X^_^  ; and 

if  x^  < X^^,  then  X^  • X^_^  - b 

The  background  ASF  estimate  for  cell  i is  X^.  If  acb,  the  median  value 
is  estimated.  Other  percentile  levels  can  be  estimated  by  adjusting  the 
ratio  of  a to  b.  The  stability  and  response  of  the  estimator  are 
controlled  by  the  magnitudes  of  a and  b.  For  nonflat  spectra,  the 
technique  is  modified  by  multiplying  by  &_  instead  of  adding,  and  by 
dividing  by  l)  instead  of  subtracting.  This  can  be  interpreted  as 
adding  and  subtracting  on  a log  scale.  This  technique  has  been 
used  extensively  at  ARL:UT  for  applications  such  as  LOFARGRAM 
prewhitening,  but  has  not  been  applied  in  this  study. 

(S)  Single  cell  detections  were  based  on  the  usual  Neyman-Pearso.. 

criteria,  wherein  the  detection  threshold  is  determined  by  the  distribu- 
tion of  the  noise  and  the  desired  false  alarm  rate.  Once  again  the 
noise  was  assumed  to  follow  a chi-square  distribution,  and  the  false 

alarm  rate  (Pe  ) was  set  to  10  J.  The  test  cell  (x. ) was  said  to  be 
ra  1 

detected  (contain  a signal)  if 


» 


where  y^  is  the  background  ASF  SPL  estimate  and  is  termed  the 
detection  coefficient.  The  detection  coefficient  is  related  to  the 
detection  threshold  (DT)  by 


(S) 


DT  - 10  log[1.5v(Cd-l)] 

where  w i.6  the  frequency  spacing  between  cells,  and  1.5  is  used  to 
normalize  the  frequency  spacing  to  its  equivalent  noise  bandwidth 
when  the  spectrum  has  been  Banned. 

(C)  The  value  of  the  detection  coefficient  is  computed  es  a function 

of  the  probability  of  false  alarm  (P-fl)  and  the  number  of  equivaler 
degrees  of  freedom  (v)  by  (Ref.  19). 

VPfa-r>  ■ {l  - £ + ZU-P , 

where  Z(P)  is  a rational  approximation  for  standard  'lorxnal  probabilities 
such  as  (Ref.  19).  This  result  is  shown  in  Fig.  1-25  for  selected  ALI 
intervals . 

(U)  The  minimum  detectable  signal  level  (MDS)  was  predicted  by  the 

methods  of  Ref.  20.  These  predictions  are  tabulated  with  the  single 
cell  processor  parameters  in  Table  1-6. 

(S)  The  false  alarm  rate  for  this  study  was  initially  specified  as 

10  , but  was  later  changed  to  10  to  provide  better  detections  against 

siioulated  third  generation  targets.  If  the  false  alarm  rate  had  been 
left  at  10  , the  standard  and  vernier  detection  thresholds  and  MDS 

would  have  been  raised  by  0.9  dB  and  1.1  da,  respectively. 

(C)  The  detection  decision  was  made  in  parallel  on  all  beaas  of  a 

frequency  cell.  If  a signal  was  detected  on  any  beacr,  the  following 
cell  parameters  were  output: 
cell  number, 

background  ASF  SPL  estimate  lor  each  beam,  and 
signal  SPL  estimate  for  each  beam. 
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TABLE  1-6 

cw  PROCESSOR  PARAMETERS  (U) 


Parameter Standard  Resolution 

ALI  Interval  5 min 

3 dB  Bandwidth  0.2197  Hz 


Number  of  Equivalent 

D<  ees  of  Freedom  (maximum) 

Noise  Window  Width 

Clipper  Probability 

Replacement  Probability 

Expected  Estimator  Bias 

False  Alarm  Probability 

Detection  Threshold 

Predicted  MDS 


(P 


det 


50%) 


161.6 

32  cells  (4.88  Hz) 
0.90 
0.50 
-0.1  dB 
10-3 

-10.6  dB/ /Hz 
-10.3  dB  (1  Hz) 


Vernier  Resolution 
5 min 

0.0176 


13.4 

64  cells  (0.78  Hz) 
0.90 
0.50 
-0.4  dB 
10~3 

-15.2  dB/ /Hz 
-14.8  dB  (1  Hz) 


(C)  where  the  signal  SPL  was  estimated  as  the  unclipped  SPL  of  the  test 
cell  minus  the  background  ASF  SPL  estimate  on  the  same  beam.  These 
results  were  stored  on  tape  for  input  to  the  line  related  processor. 

(C)  First,  the  line  related  processor  (Fig.  1-26)  sequentially  inputs 

the  detected  cell  parameters.  Based  on  these  parameters,  the  detected 
cells  for  each  ALI  interval  were  then  formed  into  lines,  and  the  line 
parameters  were  estimated.  Finally,  the  lines  were  tracked  from  each 
ALI  interval  to  the  next,  and  the  resulting  line  histories  were 
stored  on  tape.  The  navigation  reconstruction  (Ref.  21)  for  each 
ALI  interval  was  merged  with  the  line  histories  to  be  readily  available 
to  the  editing  and  display  processors.  Each  of  the  algorithms  for  line 
formation,  parameter  estimation,  and  tracking  will  now  be  described. 

(S)  Line  formation  is  the  procedure  whereby  the  detected  cells  from  an 

ALI  interval  are  grouped  together  as  lines  on  the  basis  of  the  cell 
parameters.  Ideally,  all  of  the  detected  cells  with  signals  emanating 
from  the  same  source  (i.e.,  projector  frequency,  auxiliary  equipment, 
or  propulsion  system  harmonic)  will  be  formed  into  a single  line.  The 
parameters  of  intercell  frequency  proximity  and  bearing  proximity 
(multibeam  sensors  only)  were  used  as  the  basis  for  line  formation 
decisions.  The  S/N  of  each  detected  cell  could  also  have  been  used. 
However,  the  algorithm  initially  specified  for  use  in  this  study  (Ref.  9) 
ignored  this  clue;  therefore,  all  successive  iterations  on  the  algorithm 
also  ignored  it. 

(S)  The  initial  algorithm  specified  that  contiguous  detected  cells  be 

clustered  together  as  cell  groups.  A bearing  gate  was  then  to  be  applied 
to  each  cell  group.  This  bearing  gate  would  cause  each  cell  group  to  be 
formed  into  one  or  more  lines.  The  bearing  gate  was  specified  as  *35° 
for  standard  resolution  data  and  ±65°  for  vernier  resolution  data. 

(S)  The  performance  of  the  above  algorithm  was  evaluated  using  known 

projector  lines.  It  was  found  that  a single  projector  line  frequently 
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(S)  resulted  in  multiple  formed  lines.  Further  investigation  revealed  that 
a single  projector  line  would  frequently  result  in  multiple  cell  groups 
(noncontiguous  detected  cells)  and  that  the  bearing  estimates  were  less 
stable  than  the  bearing  gates  permitted.  There  are  several  candidate 
causes  for  these  effects,  including  Doppler  smearing  during  the  5 min 
ALI,  multipath  Doppler  differences  and  frequency  interference,  source 
instabilities,  and  estimator  errors  due  to  the  low  S/N  of  some  of  the 
detected  cells  in  each  group.  The  formation  of  multiple  lines  for  each 
projector  line  resulted  in  poor  line  tracker  performance.  The  duration 
of  each  of  the.  formed  lines  was  short,  resulting  in  short  holding  times, 
even  when  continuous  detections  were  obtained. 

(S)  This  "overformation"  of  projector  lines  was  eliminated  by  broadening 

the  frequency  and  bearing  gates.  A brief  survey  of  the  data  resulted  in 
a frequency  gate  of  12  cells  for  vernier  data,  which  was  followed  by  a 
bearing  gate  whose  width  was  dependent  on  the  intercell  frequency 
spacing  (Fig.  1-27).  The  final  line  formation  algorithm  was  the  following. 

(S)  Consecutive  detected  cells  are  formed  into  a single  line  if 

(1)  the  frequency  difference  between  consecutive  detections  is 
less  than  or  equal  to  12  cells  for  vernier  data  or  two  cells  for 
standard  data,  or 

(2)  if  bearing  estimates  are  available,  the  bearing  difference  in 
degrees  between  consecutive  detections,  times  the  frequency  difference 

in  cells,  must  be  less  than  or  equal  to  90. 

The  bearing  estimates  were  computed  via  the  method  described  in 
section  II,  and  were  corrected  for  magnetic  variation  by  the  addition  of 
approximately  340°.  The  bearing  difference  was  computed  by 

AB  = min  ( | 8^  | , 360  - | B^Bj  j ) 

(C)  The  revised  line  formation  algorithm  was  found  to  perform  well  on 

projector  lines,  in  that  multiple  lines  were  seldom  formed.  This 
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FIGURE  1-27 

LINE  FORMATION  BEARING  DIFFERENCE  TOLERANCE 
AS  A FUNCTION  OF  FREQUENCY  DIFFERENCE  (C) 

AS- 79-4  7 
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(C)  reduction  in  line  loading  clearly  improved  line  tracker  performance. 

The  broadening  of  the  frequency  and  bearing  gates  also  caused  an 
increase  in  the  false  association  rate.  One  possible  effect  is  that 
closely  spaced  lines  (doublets)  may  be  formed  into  a single  line. 

(C)  An  approximation  to  the  false  association  rate  can  be  obtained 

by  modeling  the  occurrence  of  detected  cells  in  each  ALI  spectrum  as  a 

Bernoulli  process,  where  each  Bernou1li  trial  is  an  independent  single 

cell  detection  decision.  Initially  asoume  that  the  success  rate  P is 

s 

equal  to  the  probability  of  false  alarm.  Since  the  number  of  successes 
in  n trials  of  a Bernoulli  process  is  described  by  the  binomial 
probability  mass  function  (PMF),  the  probability  of  not  detecting  another 
cell  within  the  next  12  cells  can  be  expressed  as 

P = Mp°(l-P  )12-0 

ND  \0/  sv  s' 


Therefore,  the  probability  that  another  detection  will  occur  within  the 
next  12  cells  and  be  incorrectly  formed  with  the  previous  detections 
(ignoring  the  bearing  gate)  is 


12 

- 1 - (1-P  > 

s 

-3 

If  the  false  alarm  rate  were  chosen  to  be  10  , the  false  association 

rate  would  be  0.012.  This  implies  that  only  one  line  out  of  100  would 
be  incorrectly  grouped.  To  verify  this  model,  the  distribution  of 
cell  gaps  between  detected  cells  was  measured  and  fitted  with  the  appro- 
priate binomial  PMF  (Fig,  1-28).  The  corresponding  success  (detection) 
and  false  association  rates  are 
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(C) 

(A) 

P - 0.004 
s 

Pp  - 0.047 

1 

(B) 

- 0.002 

■ 0.024 

> 

(C) 

= 0.0145 

- 0.161 

, and 

(D) 

- 0.005 

- 0.058 

• 

The  false  association  rates  for  (C)  and  (D)  were  later  reduced  by 
application  of  bearing  gates,  and  the  effect  of  any  residual  error  was 
minimized  by  selecting  line  parameter  estimators  which  were  insensitive 
to  them. 

(C)  The  following  parameters  were  estimated  for  each  beam  (i)  of  a 

formed  line.  (Integrals  are  actually  summations  over  the  detected  cells.) 

• NOISE  LEVEL  = /*N^(f )dfX/*df , where  N^f)  is  the  mean 

background  ASF  level  for  each  detected  cell 
of  beam  i. 

• MAXIMUM  S/N  = maximum  S/N  of  any  detected  cell  of  beam  i = 

the  maximum  over  all  detected  cells  on  beam  i 
of  | signal  SPL/ [N^f )/ (1.5  Af)]},  where  Af  is 
the  frequency  spacing  between  Fourier 
coefficients,  and  1.5  Af  is  the  equivalent 
noise  bandwidth  of  a Manned  spectra. 

• MAXIMUM  SIGNAL  LEVEL  - signal  SPL  of  the  mrxiouo  S/N  cell  of  beam  i 

normalized  by  1.5  to  correct  for  the  Hanning 
window. 

• TOTAL  SIGNAL  LEVEL  » /s^fjdf,  where  S^f)  is  the  signal  level 

estimate  for  each  detected  cell  of  beam  i. 

• TOTAL  S/N  - TOTAL  SIGNAL  LEVEL/NOISE  LEVEL. 

The  maximum  estimates  are  used  to  address  detection  issues.  The  total 
estimates  are  similar  to  those  used  in  the  previous  ARL:UT  FVT  study 
and  to  those  predicted  by  most  models. 
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(C) 


The  following  parameters  were  estimated  for  each  formed  line. 

• POWER  MEAN  FREQUENCY  = ft*  (f  (f  )df , where  j is  the  beam 

with  the  maximum  S/N  for  the  detected  cell  of 
frequency  f. 


. rms  BANDWIDTH 


- [/(f-F)2S  (f)df// ( f ) d f ] 1/2 
= [/f2S  (f)df//S  (f)df-F2] 


where  F is  the  POWER  MEAN  FREQUENCY. 


• SIGNAL  BEARING 


arctan 


/*sin[B  (f)]*S 
/coslBs(f))xSj(f)df  \ 


where  B (f)  is  the  noise  bearing  estimate  of  the  detected  cell 
s 

at  frequency  f. 


• NOISE  BEARING 


(/sin(Bn(f))/si(f)dfi 

|ycos[Bn(f) 


where  Bn(f)  is  the  noise  bearing  estimate  of  the  detected  cell 
at  frequency  f. 

Note  that  the  line  bearing  estimators  are  simply  weighted  summations  of 
the  cell  bearing  estimates,  where  the  weight  is  the  signal  SPL  on  the 
beam  with  maximum  S/N,  or  its  inverse.  The  detected  cell  noise  bearing 
estimates  are  computed  in  the  same  manner  as  the  signal  bearing  estimates. 


(C)  Any  received  signal  is  expected  Co  occupy  more  than  one  frequency  cell. 
Even  if  the  analysis  bandwidth  were  made  sufficiently  coarse  and  the  AH 
interval  sufficiently  short  so  that  Doppler  shifts,  multipath  effects,  and 
source  instabilities  were  contained  in  a band  whose  width  was  less  than 
the  analysis  resolution,  the  received  signal  would  still  occupy  more  than 
one  cell  because  of  the  sidelobo  leakage  of  the  FFT.  This  suggests  that 
signal  detection  should  be  performed  by  techniques  which  utilize  this 
phenomenon.  Consider  for  example  a detector  where  two  adjacent  cells 
arc  required  to  exceed  a threshold.  The  threshold  would  be  determined 
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(C)  so  as  to  maintain  the  desired  false  alarm  rate  against  stationary 
Gaussian  input  data.  This  determination  would  need  to  account  for 
the  cross frequency  correlation  induced  by  overlapped  FFT  processing. 
The  resultant  threshold  would  be  lower  than  a single  cell  threshold 
with  the  same  false  alarm  rate.  This  technique  could  be  generalized 
into  a matched  bandwidth  detector  which  would  contain  multiple 
bandwldths  and  their  associated  thresholds.  The  MDS  of  the  matched 
bandwidth  detector  would  be  significantly  lower  than  that  of  a single 
cell  detector.  Such  a technique,  for  detecting  lines  instead  of 
single  cells,  would  also  eliminate  the  necessity  for  line  formation, 
and  should  lead  to  better  line  parameter  estimates. 

(C)  The  final  function  of  the  line  related  processor  (Fig.  1-26)  was 
line  tracking.  The  estimated  line  parameters  were  used  to  link  the 
lines  detected  during  the  current  ALI  interval  to  the  current  set  of 
active  line  histories.  If  the  line  could  not  be  linked,  a.  new  line 
history  was  created.  Line  histories  weve  purged  from  the  active  set 
after  failing  to  be  linked  during  a specified  number  of  contiguous 
ALI  intervals. 

(S)  The  initial  line  tracker  specification  (Ref.  9)  used  frequency  and 
bearing  proximity  as  the  basis  for  a linking  decision.  For  linking, 
the  current  line  was  required  to  be  within  a frequency  gate  and  a 
bearing  gate  (i 30°)  of  an  active  line  history.  The  frequency  gate  was 

0.2  Hz  0.01  * (center  frequency  of  line  history) 

for  standard  resolution,  and 

+ 1 cell,  - 2 eells  in  the  50  Hz  region, 

+ 2 cells,  - 4 cells  in  the  150  Hz  region,  and 

+ 2 cells,  - 8 cells  in  the  300  Hz  region 

for  vernier  resolution.  After  a line  was  linked,  the  line  history 

frequency  and  bearing  were  to  be  updated  by 

74 


(S) 


updated  frequency  ■ (old  frequency  + new  frequency) /2 


and 

updated  bearing  ■ (old  bearing  + new  bearing) /2. 

(C)  When  applied  to  known  cw  signals,  the  performance  of  this  line 
tracker  was  found  to  be  poor.  Even  though  the  projectors  were  being 
towed  at  low  speed  (5  kt),  the  line  tracker  gates  were  insufficiently 
wide  to  link  their  detections.  The  result  was  that  several  short  line 
histories  were  formed  for  each  projector  line.  Line  tracker  performance 
measures  (i.e.,  holding  time)  were  not  evaluated  by  this  study,  but  a 
fair  performance  level  was  required  so  that  the  desired  measures  (i.e., 
percentage  detection)  could  be  obtained.  To  provide  this  performance 
level,  a new  tracking  algorithm  was  evolved.  The  sole  selection 
criteria  applied  to  this  algorithm  was  its  performance  against  known 
projector  lines. 

(C)  The  revised  line  tracking  algorithm  consists  of  a distance  function 
and  a four-step  link  selection  procedure.  The  distance  function  was 
designed  to  provide  a measure  of  the  similarity  between  a line  and  a 
line  history.  The  parameters  used  by  this  measure  were  frequency 
proximity,  bearing  proximity  (if  available),  and  S/N.  Other  clues  such 
as  bandwidth  and  line  dynamics  were  not  investigated.  A general  purpose 
link  selection  procedure  was  first  devised,  and  then  used  to  evaluate 
several  distance  functions. 

(C)  During  the  first  step  of  the  link  selection  procedure,  the  distance 

between  each  formed  line  and  each  line  history  was  computed.  If  the 

distance  was  less  than  some  fixed  maximum  (D  ),  the  line  history  was 

max 

added  to  the  list  of  candidate  links  for  that  formed  line.  As  many  as 
seven  candidate  links  were  allowed  for  each  line.  During  the  second 
step,  the  preferred  link  was  selected  from  the  list  of  candidate  links 
for  each  sensor.  The  highest  preference  was  given  to  the  line  history 
with  the  most  detections  (longest).  If  the  longest  histories  were  of 
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(C)  equal  length,  the  closest  (smallest  distance)  history  was  preferred. 

If  the  histories  were  equally  close,  the  highest  frequency  history  was 
preferred.  This  ordering  was  evolved  through  many  iterations  on  the 
data.  The  longest  history  criterion  is  especially  essential  for  the 
proper  tracking  of  known  cw  signals.  For  example,  because  of  a course 
change  within  an  ALI  interval,  two  lines  were  sometimes  formed  at 
significantly  different  frequencies  from  a single  projector  line.  If 
the  longest  history  criterion  was  not  applied  during  the  following 
ALI  interval,  the  precourse  change  history  would  not  be  linked  to 
postcourse  change  detections.  Two  examples  of  this,  during  a single 
projector  tow,  will  be  shown  in  Fig.  1-31. 

C)  During  the  third  step  of  the  link  selection  procedure,  conflicts 
are  resolved  between  lines  having  the  same  preferred  line  history. 
Preference  is  given  to  the  closest  line.  Finally,  when  all  of  the 
conflicts  have  been  resolved,  the  fourth  step  is  performed.  All  of 
the  formed  lines  for  which  links  have  been  selectee  are  now  associated 
with  the  specified  line  histories.  The  line  history  parameters 
(frequency,  baLdwidth,  S/N,  and  bearing)  are  updated  by 

updated  value  * (old  value  + new  value) /2  , 


Those  formed  lines  for  which  no  link  was  selected  are  used  to  create 
new  line  histories. 


Line  histories  were  purged  from  the  active  list  after  six  consecutive 
ALI  intervals  during  which  no  new  links  were  generated.  Thu  active  list 
had  a capacity  of  400  line  histories. 

The  form  of  the  distance  function  was  designed  to  meet  the  anticipated 
worst  case  tracker  requirements  for  these  data.  When  the  projector 
platform  was  proceeding  along  a straight  line  radial  track  at  a significant 
range  (>10  nei)  from  the  receiver,  the  frequency  and  bearing  gates  were 
required  to  be  only  as  wide  as  the  frequency  and  bearing  estimator  errors. 
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(C)  Since  the  track  was  usually  nonradial,  the  gates  ver*  widened  to  account 
for  Doppler  and  bearing  rates.  Because  of  the  low  platform  speed 
(<5  kt),  these  rates  were  small  and  could  have  been  tracked  by  the 
initial  algorithm. 

(C)  The  tow  platforms  underwent  frequent  course  changes,  sometimes  in 
excess  of  120°  (see  Volume  IV).  Such  course  changes  caused  very  large 
and  rapid  Doppler  frequency  shifts  while  leaving  the  bearing  relatively 
constant.  The  normalized  frequency  shift  for  a course  reversal  along 
a 5 kt  radial  tow  is  0.34%  of  the  signal  frequency,  where  the  effect 
of  receiver  depth  has  been  ignored.  The  revised  line  tracker  was 
therefore  provided  with  a broad  frequency  gate  whose  width  is  proportional 
to  the  line  history  center  frequency. 

(C)  The  tow  platforms  frequently  passed  within  close  proximity  of  the 

receiver  sites.  Unlike  course  changes,  both  the  Doppler  frequency  rate 

and  bearing  rate  can  become  large  at  CPA.  The  maximum  Doppler  frequency 

rate  occurs  when  the  platform  is  directly  above  the  receiver,  and  is 
2 

equal  to  V /rc  where  V is  the  platform  speed,  r is  the  range  at  CPA, 
and  c is  the  velocity  of  sound.  For  a 3 kt  speed,  1525  m (5000  ft)  deep 
receiver,  and  a 5 min  AH,  the  Doppler  shift  becomes  approximately  0.09%  of 
the  signal  frequency.  The  bearing  estimate  could  change  by  180*  between 
two  consecutive  ALI  intervals  about  CPA.  In  addition,  the  bearing 
estimates  become  less  stable  when  the  platform  is  directly  above  the 
D1FAR  sensor. 

(C)  The  large  frequency  shifts  and  small  bearing  shifts  due  to  course 
changes,  and  the  smaller  frequency  shifts  and  large  bearing  shifts  at 
CPA,  suggested  that  the  line  cracker  use  a distance  function  of  the 
following  form. 


D * a iii/i  * &Ab 
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(C)  where  Af  and  Ab  are  the  frequency  (hertz)  and  bearing  (degree) 

differences  between  a line  and  its  candidate  line  history,  f is  the 
frequency  (hertz)  of  the  line  history,  and  a and  3 are  constants.  This 
function  allows  one  gate  to  broaden  as  the  other  narrows,  thus  matching 
the  anticipated  conditions, 

(C)  When  values  ..or  ct  ana  3 were  chosen  to  permit  the  line  tracker  to 
perform  well  with  the  known  projector  lines,  the  performance  against 
clutter  wfcs  extremely  poor.  Unrelated,  spurious  detections  were 
tracked  as  lines,  resulting  in  increased  processor  loading.  The 
predominant  characteristic  of  these  false  line  histories  was  that  they 
were  composed  of  lines  with  very  low  S/N.  The  form  of  the  distance 
function  was  therefore  modified  to  discriminate  against  such  links  by 
narrowing  the  frequency  r t aring  gates  for  low  S/N  lines. 

D = oAf/f  + BAb  - y I (S.  K i 

where  S/N^  and  S/N^  are  the  signal-to-noise  ratios  (dB//Hs)  of  the  line 
and  the  candidate  line  history. 

(S)  After  several  iterations,  it  was  heuristically  determined  that  a 
frequency  difference  of  1.5  vernier  cells  at  50  Hz  was  equivalent  to 
1 dB  S/N,  which  was  equivalent  to  0,1°  bearing  difference.  The 
generalized  form  of  the  distance  function  was 

D * 6144  A^lr!  + 0.1  Ab  (deg)  - [ S/N.  (dB//Hz)  + S/N  (dB//Hz)] 

I V.HZ ) L n 

The  maximum  distances  for  which  linking  was  allowed  (D  ) were  chosen 

max 

as  36,  if  bearing  was  available,  and  30,  if  not.  For  example,  lines  on 
single  beam  (no  bearing)  sensors  with  vernier  threshold  levels  (-15  dB) 
were  not  allowed  any  frequency  difference,  whereas  high  level  signals 
(0  dB)  were  allowed  frequency  rates  of  up  to  0.49%.  Little  weight  was 
placed  oii  bearing  because  of  the  rapid  changes  expected  at  short  range. 
The  value  of  this  clue  would  be  increased  if  it  were  known  when  the 


(S)  source  was  at  short  range.  It  should  be  possible  to  obtain  such  a 
clue  from  the  DIFAR  estimate  of  vertical  arrival  angle  achievable 
with  the  vertical  sensor  discussed  in  section  II.  Time  did  not 
permit  this  investigation  to  be  included  in  this  study. 

(C)  The  form  of  the  line  tracker  distance  function  involves  the 

implicit  assumption  that  the  signal  should  occur  at  the  same  frequency 
and  bearing  during  each  ALI  interval.  Any  deviation  from  the  historic 
values  causes  an  increase  in  the  distance  measure.  Also,  the  historic 
values  are  updated  in  a manner  that  implies  that  only  a portion  of 
the  measured  change  is  real.  The  effect,  when  the  line  is  undergoing 
real  frequency  and  bearing  shifts,  is  to  successively  increase  the 
distance  along  the  real  link,  and  thus  increase  the  probability  of 
error.  Another  technique  that  merits  consideration  is  that  used  in  the 
BEARTRAP  program  (Ref.  22).  The  frequer  ,/  and/or  bearing  rates  could 
be  estimated  from  recent  line  history  values  and  used  to  predict  the 
future  values  of  these  parameters.  The  distance  function  would  then 
be  based  on  the  predicted  values.  In  addition  to  decreasing  the  chance 
of  error,  this  technique  furnishes  estimates  of  line  dynamics  which  are 
of  value  for  other  functions  such  as  signature  formation.  Also,  if  a 
target  was  being  tracked,  but  insufficient  lines  were  available  for 
classification,  these  parameter  estimates  could  be  used  to  key  a 
dynamic  line  integrator  (DLI)  process  (Ref.  23)  and  thereby  increase 
the  probability  of  detection  on  other  related  target  lines. 

(C)  Tne  line  tracker  distance  function  was  evaluated  for  several 

edited  projector  line  histories,  and  the  values  plotted  in  Figs.  1-29 
and  1-30.  The  lines  on  each  figure  denote  the  restrictions  on  linking 
imposed  by  the  current  function.  These  figures  reveal  that  most,  but 
not  all,  of  the  links  could  survive  more  stringent  restrictions.  The 
advantage  of  tighter  linking  criteria  is  a lower  clutter  (line/processor 
loading)  rate.  The  figures  also  reveal  that  a few  of  the  correct  links 
were  rejected  by  the  current  criteria.  These  links  were  added  in  a later 
manual  editing  stage.  Automatic  generation  of  accurate  line  histories 
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FIGURE  1-29 

INE  TRACKER  DISTANCE  FUNCTION  VALUES  FOR  THE  17  NOV  PROJECTOR  LINES 
OBSERVED  WITH  THE  SINGLE  CARDIOIDS  SENSORS  AND  VERNIER  RESOLUTION  d 
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(C)  may  also  require  a post-line-tracker  routine  which  ; ^ociates  terminated 
line  histories  with  new  histories  that  have  similar  characteristics. 

(S)  Finally,  the  line  formation  and  tracking  algorithms  that  were 

developed  for  these  data  might  not  perform  well  on  real  data.  Further 
revisions  would  be  required  to  track  high  speed  (>10  kt)  targets  and/or 
those  with  unstable  lines  (Ref.  24). 

(C)  The  line  histories  output  by  the  line  related  processor  were  manually 
edited  prior  to  the  generation  of  final  data  products.  This  editing  was 
accomplished  by  means  of  high  resolution  detection  history  plots  such  as 
Fig.  1-31.  All  of  the  detected  lines  in  a 2 Hz  band  centered  about  each 
projector  frequency  were  plotted  as  a function  of  time.  Detections 
based  only  on  a single  cell  of  a single  beam  are  denoted  by  X's. 

Multicell  and/or  multibeam  detections  are  denoted  by  0's.  Linking  is 
denoted  by  dashed  lines.  Bearing  is  denoted  by  a solid  line  emanating 
from  each  detection.  The  length  of  this  line  is  proportional  to  S/N. 

Ground  truth  range  nautical  miles  (nmi)  and  bearing  (degrees)  are  plotted 
across  the  top.  These  kinds  of  displays  were  used  to  select  which  history, 
histories,  or  segments  of  histories  corresponded  to  each  projector  line. 
Following  history  selection,  the  measured  cw  signal  parameters  were  plotted 
as  a function  of  time  in  the  format  of  Fig.  1-32.  These  line  history  plots 
and  a description  of  their  contents  can  be  found  in  Volumes  II  and  III. 
Following  this  editing  and  parameter  extraction,  detailed  and  summary  cw 
data  products  were  generated.  These  products  and  their  descriptions 
are  found  in  Volumes  II  through  IV, 

(C)  The  last  data  products  to  be  generated  by  the  cw  measurement  system, 
clutter  statistics,  were  actually  a measure  of  the  ambient  sound  field. 

As  such,  it  was  desirable  to  eliminate  any  exercise  induced  artifacts. 

For  example,  Fig.  1-33  contains  all  of  the  line  histories  from  the  single 
cardioids  sensor  at  site  A1  during  the  17  November  field  event.  Figure  1-34 
is  an  example  of  the  corresponding  clutter  display  for  these  data.  Numerous 
lines  were  held  throughout  the  12  h interval.  Closer  inspection  reveals 
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FIGURE  1-33 


^S5-FVT  DETECTION  OVERVIEW  FOR  SINGLE  CARDIOIDS  SENSC 
AT  oITE  A1  INCLUDING  ALL  EXERCISE  RELATED  LINES  (",) 
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FIGURE  1-34 

CLUTTER  TIME  SERIES  FOR  SINGLE  CFIRDIOIDS  SENSOR 
E fil  INCLUDING  FILL  EXERCISE  RELATED  LINES  (U) 
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(C)  that  many  of  these  lines  (Fig.  1-35)  emanated  from  a single  projector 
platform  (CFAV  KAPUSKACING) , even  though  only  three  lines  appear  in  the 
projector  operations  logs.  Many  of  the  remaining  lines  (Fig.  1-36) 
can  also  be  attributed  to  exercise  related  sources  (ACODAC,  CFAV  QUEST, 
and  R/V  CHAIN).  When  all  of  these  known  lines  were  eliminated  (Fig.  1-37), 
the  remaining  detections  provided  more  accurate  clutter  measurements 
(Fig.  1-38).  All  of  the  clutter  results  in  this  study  were  derived  with 
a similar  editing  procedure.  Even  after  editing,  some  exercise  artifacts 
remain  to  bias  the  clutter  measurements.  The  many  high  frequency  lines 
of  Fig.  1-37  beginning  at  321/2230  were  due  to  the  GNATS  system  aboard 
USS  GREENLING  (SSN  614).  Additional  GNATS  lines  can  be  observed  at 
321/1300  as  USS  GREENLING  passed  CPA.  However,  even  if  some  residual 
artifacts  remain,  the  edited  clutter  measurements  are  much  more  repre- 
sentative of  the  true  values  than  the  unedited  measurements. 
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HSS-FVT  DETECTION  OVERVIEW  TOR  SINGLE  CARDIOIDS  SENSOR 
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VII.  CONCLUSIONS  AND  RECOMMENDATIONS 


(C)  A measurement  system,  typical  of  those  envisioned  for  use  by  MSS, 
has  been  developed  at  ARL:UT.  The  design  was  based  on  the  current 
ARL:UT  measurement  system  and  on  the  guidelines  provided  by  the  sponsor. 
The  algorithms  were  finalised  after  many  trials  and  were  evaluated  on 
the  basis  of  their  performance  against  known  projector  lines.  While 
clearly  suboptimal,  the  performance  of  these  algorithms  was  adequate 
for  the  purposes  of  this  study.  Several  improvements  to  these  algorithms 
were  recommended  for  further  study,  including  the  following. 

• Fonnation  of  a sensor  with  a vertical  dipole  response  from  a 
single  DIFAR  array. 

• A recursive  technique  for  background  ASF  estimation  with  improved 
accuracy  and  efficiency. 

• A matched  bandwidth  detector  technique  with  a lower  MDS  than  that 

of  the  single  cell  detector  currently  in  use. 

, • • ** 

• Use  of  a linear  predictor  by  the  line  tracker  to  provide  more 
accurate  tracks  and  estimates  of  line  dynamics  for  signature 
formation. 

• A technique  for  vertical  arrival  angle  estimation  from  a single 
DIFAR  array,  which  in  turn  would  furnish  a range  estimate  to 
improve  line  tracker  and  signature  formation  performance . 

(U)  Finally,  future  sensor  performance  analysis  will  not  necessarily 
require  this  type  of  processor,  except  to  obtain  clutter  results. 

Rather,  a cw  measurement  technique  which  takes  advantage  of  the  known 
source  and  source  platform  parameters  would  more  readily  provide 
projector  line  histories.  Such  histories  would  include  measurements  at 
low  S/N.  many  of  which  would  not  be  detected  by  the  current  processor. 
These  histories  could  then  be  used  to  predict  what  would  have  been 
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(U)  observed  by  other  types  of  processors.  These  results  would  be  more 
complete  and  more  easily  obtained, 
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APPENDIX  I-A 

STANDARD  RESOLUTION  INTERCHANNEL  PHASE  DIFFERENCES 
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(U)  The  following  set  of  plots  contains  the  measured  interchannel  phase 

differences  as  measured  at  the  output  of  the  ARL:UT  measurement  system. 

The  first  two  plots  contain  measurements  of  the  ;'RL:UT  measurement  system 
alone.  The  first  set  of  measurements  was  performed  just  prior  to  digiti- 
zation of  the  analog  data;  the  second  set  was  performed  immediately 
following  this  digitization.  The  remaining  plots  contain  before  and  after 
measurements  based  on  the  header  calibration  signals  of  the  government 
furnished  (GFE)  analog  tapes. 
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APPENDIX  1-B 

STATISTICAL  STABILITY  OF  AMBIENT  SOUND  FIELD  MEASUREMENTS 
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(U)  An  estimate  of  the  statistical  stability  cf  the  ambient  sound  field 

(ASF)  measurements  reported  in  this  study  can  be  obtained  by  modeling 
the  ASF  as  a stationary  Gaussian  random  process.  The  out’ut  of  a FFT 
processor,  such  as  that  used  herein,  for  such  input  will  follow  a distri- 
bution which  can  be  closely  approximated  by  a chi-square  distribution 
(Ref.  15).  The  chi-square  distribution  is  completely  specified  by  a 
single  parameter,  the  number  of  degrees  of  freedom.  In  this  appendix, 
results  will  be  shown  regarding  the  number  of  equivalent  degrees  of 
freedom  possessed  by  an  estimate  of  the  average  amount  of  power  of  a 
Gaussian  random  noise  process  in  a given  frequency  band,  where  the  estimate 
is  to  be  obtained  by  combining  adjacent  bins  of  a higher  resolution  power 
spectrum  estimate.  (The  number  of  equivalent  degrees  of  freedom  (EDF) 
of  an  estimate  refers,  as  usual,  to  tho  chi-square  approximation  to  the 
actual  distribution  of  the  estimate.)  The  results  are  derived  in 
Ref.  16,  where  it  is  shown  that,  when  the  high  resolution  spectral 
estimates  are  obtained  by  averaging  50Z  overlapped  FFTs,  the  number  of 
equivalent  degrees  of  freedom  of  an  estimate  of  the  power  in  a band 
containing  K+l  FFT  bins  is  given  by 
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(U)  where 

(1) 


/CD 

|w(f)|2df  , 


(2)  Ij(M)  = |/  W(f)W*(f^k>£)df 


(3)  I2(k,i) 


/•’ 

/ W(f)W 


(f+Ak^)ejirLf  df 


(4)  P is  the  number  of  FFTs  used  to  form  the  high  resolution  spectral 
estimate, 

(5)  L is  the  length  of  each  FFT  in  seconds. 


(6) 


k-£ 


M L * 

(7)  a k aK  are  we*Shting  coefficients  assigned  to  the  FFT 

' 2 2 

(8)  bins  which  are  summed  to  produce  the  estimate,  and 

(9)  W(f)  is  the  spectral  window  employed. 
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(U)  For  Hanned  spectra, 


wm  = 1 sinc(Lf) 
,0  zh-CLf)2] 


and  a numerical  evaluation  of  the  integrals  and  I20 
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A straightforward  computation  for  I gives 


I = 1. 406^S  x 10'1  L2 
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(U)  Substituting  back  into  Eq.  (B2),  we  have,  for  Banned  spectra, 
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where 

Co  = 1.40625  * 10' 1 ♦ (3.90625  * 10'3)  , 

Cj  = 6.25000  x 10*2  ♦ (2.81448  * 10~3) 

C,  « 3.90625  * 10'J  ♦ (9.76562  * 10'4)  . and 
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(U)  This  result  was  used  to  cooputc  the  EOF  for  each  frequency  band 

froes  which  ASF  aeasurettents  were  obtained.  Selected  percentile  levels 
denoting  the  spread  of  the  corresponding  chi-square  distribution  about 
its  scan  vaiue  were  plotted  in  Fig.  1-19. 

(U)  Equation  B(2)  was  also  used  to  predict  the  statistical  variation 

in  noise  gain  which  could  be  expected  In  the  absence  of  a signal  when 
the  noise  process  was  isotropic,  stationary,  and  Gaussian.  Here  the 
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(U)  noise  gain  Z(f)  is  defined  as  the  ratio  x(f)/Y(f),  where  x(f)  is  an 
estimate  of  the  noise  power  contained  in  a frequency  band  about  f as 
measured  by  a directional  receiver,  and  Y(f)  is  the  corresponding 
estimate  for  a nondirectional  receiver.  In  this  case,  50%  overlapped 
Banned  FFTs  were  used  to  form  estimates  of  the  power  spectra  for  both 
the  directional  and  nondirectional  receivers.  It  was  desired  to 
convert  the  standard  uniform  resolution  of  these  spectral  estimates 
to  logarithmic  resolution  by  dividing  the  frequency  range  into 
1/10-octave  bands  and  summing  all  the  FFT  bins  in  each  band  (with  uniform 
weighting  across  the  band).  Using  Eq.  (B2),  the  number  of  degrees  of 
freedom  v for  the  approximating  chi-square  distribution  could  be 
determined  for  both  x(f)  and  Y(f)  for  each  1/10-octave  band.  Then, 
assuming  the  directional  and  nondirectional  measurements  were  statis- 
tically independent,  standard  techniques  can  be  applied  to  show 


■i»«  • ■)  • ST 


f (x)f (v)dxdy 


o o 


where 


x > 0 

X % 0 


This  integral  was  evaluated  numerically  for  several  values  of  s,  and 
the  results  were  also  plotted  in  Fig.  1-19. 
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APPENDIX  I-C 

BIAS  OF  BACKGROUND  AMBIENT  SOUND  FIELD  MEAN  LEVEL  ESTIMATOR 
IN  THE  ABSENCE  OF  ANY  SIGNALS 
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(C)  The  background  ambient  sound  field  (ASF)  mean  level  estimator  used 

by  the  cw  signal  processor  in  this  study  was  based  on  a clipper/ 
replacement  technique.  The  sound  pressure  level  (SPL)  of  each  frequency 
cell  which  exceeds  a clipping  threshold  is  replaced  by  an  estimate  of 
the  medium  SPL.  A frequency  band  of  the  clipped  spectrum  is  then  averaged 
to  estimate  the  mean  background  ASF  SPL.  When  no  signals  are  present, 
the  clipping  of  high  level  ASF  levels  will  result  in  the  background 
estimator  being  biased  low.  The  amount  of  bias,  for  stationary  Gaussian 
input  data,  will  now  be  derived. 

(U)  When  stationary  Gaussian  data  is  processed  by  a FFT  technique  such 

as  that  used  in  this  study,  the  distribution  of  the  output  spectra  is 

2 

closely  approximated  by  the  product  of  a chi-square  distribution  (x  ) 
and  a scaling  term  which  can  be  neglected  in  what  follows.  The  mean 
value  of  the  distribution  is 


E[x(£) ] 


(x)dx 


* Y 


where  y is  the  number  of  degrees  of  freedom  defining  the  chi-square 
distribution. 
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Define  the  clipping  threshold  (T  ) and  replacement  value  T„  as 

C K 


and 


tr  = yCr 


where  Cc  and  CR  are  termed  the  clipper  and  replacement  coefficients  and 

are  chosen  to  set  the  clipper  threshold  and  replacement  values  at  the 

desired  probabilities  (P  and  P ).  The  expected  value  of  the  clipped 

c R 

spectrum  is 


E(xc(f)] 


/ C x X^(x)dx  * f Tr  X^(x)dx  , 


where 


X^(x) 


[‘Y/2r(i-)]'1 


xV/2-le-x/2 


The  bias  of  the  ASF  estimate  is 


bias  * 


Efxc(f)] 

TRfTT 


[y2Y/2r(v)]  J C XY/2e“x/idx  + CR(1-Pc) 


For  values  of  v which  are  even  integers,  the  bias  bocoiaos 


126 

CONFIDENTIAL 


CONFIDENTIAL 


(C)  bias 
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(y/2) * T K 
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(Y/2)! 
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+ CK<1'Pc) 


(U)  This  function  is  plotted  in  Fig.  1-24  as  a function  of  Y and  for  a 
replacement  coefficient  corresponding  to  a probability  of  50%. 
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9 OP-951F 

1°  OP-951FI 

11  OP-981 

Commander 

Naval  Air  Development  Center 
Department  of  the  Navy 
Warminster,  PA  18974 

12  Attn:  Code  3091 

13  Code  303 

Commander 

Naval  Air  Systems  Command 
Department  of  the  Navy 
Washington,  DC  20360 
Attn:  Code  PMA-264A 
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15  Commanding  Officer 

Naval  Intelligence  Support  Center 
4301  Suitland  Road 
Washington,  DC  20390 
Attn:  Code  222 

16  Project  Manager 

Anti-Submarine  Warfare  Systems  Project  Office 
Department  of  the  Navy 
Washington,  DC  20360 
Attn:  ASW  01T 

Commanding  Officer 

Naval  Ocean  Research  and  Development  Activity 
Department  of  the  Navy 
NSTL  Station,  MS  39529 

17  Attn:  S.  Marshall,  Code  340 

18  G.  Lewis,  Code  500 

19  - 20  Commanding  Officer  and  Director 

Defense  Documentation  Center 
Cameron  Station,  Building  5 
5010  Duke  Street 
Alexandria,  VA  22314 

21  TRW,  INC. 

TRW  Defense  and  Space  Systems  Group 

Washington  Operation 

7600  Colshire  Drive 

McLean,  VA  22101 

Attn:  W.  Morley 

22  Sanders  Associates,  Inc. 

95  Canal  Street 
Nashua,  NH  03606 

Attn:  L.  Gagne 

23  Office  of  Naval  Research 
Resident  Representative 
Room  508,  Federal  Building 
Austin,  TX  78701 

24  Glen  £,  Ellis,  ARL : UT 

25  Loyd  D,  Hampton,  ARL:UT 
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26 

Kenneth  E.  Hawker,  ARL:UT 

27 

Stephen  K.  Mitchell,  ARLsUT 

28 

Clark  S.  Penrod,  ARL:UT 

29 

Jack  A.  Shooter,  ARL:UT 

30 

Library,  ARL:UT 
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DEPARTMENT  OF  THE  NAVY 

OFFICE  OF  NAVAL  RESEARCH 
875  NORTH  RANDOLPH  STREET 
SUITE  1425 

ARLINGTON  VA  22203-1995 

IN  REPLY  REFER  TO: 


5510/1 

Ser  32 1OA/0 11/06 
31  Jan  06 


MEMORANDUM  FOR  DISTRIBUTION  LIST 

Subj:  DECLASSIFICATION  OF  LONG  RANGE  ACOUSTIC  PROPAGATION  PROJECT 
(LRAPP)  DOCUMENTS 

Ref:  (a)  SECNAVINST  5510.36 

Enel:  (1 ) List  of  DECLASSIFIED  LRAPP  Documents 

1 . In  accordance  with  reference  (a),  a declassification  review  has  been  conducted  on  a 
number  of  classified  LRAPP  documents. 

2.  The  LRAPP  documents  listed  in  enclosure  (1)  have  been  downgraded  to 
UNCLASSIFIED  and  have  been  approved  for  public  release.  These  documents  should 
be  remarked  as  follows: 

Classification  changed  to  UNCLASSIFIED  by  authority  of  the  Chief  of  Naval 
Operations  (N772)  letter  N772A/6U875630,  20  January  2006. 

DISTRIBUTION  STATEMENT  A:  Approved  for  Public  Release;  Distribution  is 
unlimited. 

3.  Questions  may  be  directed  to  the  undersigned  on  (703)  696-4619,  DSN  426-4619. 


^ — Jr. 

BRIAN  LINK 
By  direction 
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PEO  LMW  Det  San  Diego  (PMS  181) 

DTIC-OCQ  (Larry  Downing) 

ARL,  U of  Texas 

Blue  Sea  Corporation  (Dr.Roy  Gaul) 

ONR  32B  (CAPT  Paul  Stewart) 

ONR  3210A  (Dr.  Ellen  Livingston) 

APL,  U of  Washington 

APL,  Johns  Hopkins  University 

ARL,  Penn  State  University 

MPL  of  Scripps  Institution  of  Oceanography 

WHO! 

NAVSEA 

NAVAIR 

NUWC 

SAIC 
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